m m.^=— 

QQC  FILE  COPY  AJ)A0  4 


Technical  Note 


Waveform  and  Clutter  Filter  Design 
for  a Satellite-Borne  Bistatic  Radar 
for  Aircraft  Detection 


I’rcrnred  for  the  Oefcnse  Advanced  Research  Projects  Agency 
under  Plectrmic  Systems  Division  Contract  F19628-76»C*0002  by 

Lincoln  Laboratory 

MASSACHUSETTS  INSTITUTE  OF  TECHNOLOGY 
LuxiNtiTON,  Massachusetts 


Approved  for  public  releisei  distribution  unlimited. 


\ ) 


! 


The  work  reported  in  this  document  was  performed  at  Lincoln  Laboratory, 
a center  for  research  operated  by  Massachusetts  Institute  of  Technology. 
This  work  was  sponsored  by  the  Defense  Advanced  Research  Projects 
Agency  under  Air  Force  Contract  F19628-76-C-0002  (ARPA  Order  600). 

This  report  may  be  reproduced  to  satisfy  needs  of  U.S.  Government  agencies. 


The  view  and  conclusions  containr-d  in  this  document  are  those  of  the 
conL-ctor  and  should  net  be  interpreted  as  necessarily  rep-esenting  the 
official  policies,  either  expressed  or  implied,  of  the  United  States 
Government. 


llsis  lechnica!  leport  has  been  reviewed  and  is  approved  for  publication. 
FOR  THE  COMMANDER 


Rg^ond  L.  L^jtselie,  Lt.Co!.,  USAF 
Chief,  ESD  Lincoln  Laboratory  Picject  Office 


MASSACHUSETTS  INSTITUTE  OF  TECHNOLOGY 


LINCOLN  LABORATORY 


WAVEFORM  AND  CLUTTER  FILTER  DESIGN 
FOR  A SATELUTE-BORNE  BISTATIC  RADAR 
FOR  AIRCRAFT  DETECTION 


R.  A.  GOVNDRY 
Group  37 

C.  J.  KVRYS 
Group  31 


TECHNICAL  NOTE  1977-7 
28  JUNE  1977 


Approved  for  public  release;  distribution  unlimited. 


D D C 


B 


LEXINGTON 


MASSACHUSETTS 


I 


p 

I 


j 

1 


I 


i 


ABSTRACT 


Waveform  and  clutter  filtering  for  the  detection  of 
moving  aircraft  in  the  presence  of  ground  clutter  by  c space- 
borne  bistatic  radar  are  described.  The  waveform  deblc^ns  are 
based  on  mappings  of  the  aircraft  and  clutter  onto  in.j  bistatic 
range-Doppler  plane.  Low  sidelobe  digital  filter  d-'.?igns  for 
the  required  clutter  rejection,  coverage  character  ,:.tics  of  a 
bistatic  radar  on  representative  satellite  orbits,  and  system 
requirements  for  a North  American  surveillance  ff-ioo  for  bomber 
detection  are  included. 
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I. 


INTRODUCTION 


A satellite-borne  radar  system  for  aircraft  detection 
will  be  viewing  targets  against  the  background  of  the  earth 
(ground  clutter)  for  most  viewing  geometries.  A bistatic  radar 
system  may  offer  the  potential  for  target  cross  section  enhance- 
ment and  clutter  cross  section  reduction  by  virtue  of  the  viewing 
geometry.  Whether  this  potential  exists  would  have  to  be  estab- 
lished by  measurements.  A performance  measure  of  a radar  system 
is  its  ability  to  detect  targets  in  a cluttered  environment  by 
either  avoiding  or  suppressing  the  clutter.  Clutcer  suppression 
can  be  achieved  through  a judicious  choice  of  waveform  and  filter 
design.  This  report  presents  the  results  of  an  investigation  of 
waveform  design  and  clutter  filtering  techniques  for  use  in  a 
satellite-borne  bistatic  radar  for  aircraft  detection. 

System  requirements  for  detection  of  a Soviet  bomber 
attack  on  the  United  States  are  presented  and  the  potential 
advantages  of  a bistatic  radar  system  over  a monostatic  radar 
system  are  reviewed.  Estimates  of  the  power-aperture  required 
of  a satellite-borne  bistatic  radar  system  to  provide  various 
coverage  rates  are  made  to  indicate  the  relative  size  of  the 
system.  Any  satellite-borne  sensor  system  requires  communication 
and  ground  control  centers  which  make  use  of  the  information 
gathered  by  the  system.  These  issues  are  recognized  but  were 
not  investigated  since  they  do  not  impact  the  investigation  of 
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waveforms  and  clutter  filtering. 


A satellite-borne  bistatic  radar  system  requires  two 
satellites,  one  for  the  transmitter  and  one  for  the  receiver. 

Two  satellites  in  the  same  or  different  orbits  can  provide  unique 
coverage  of  certain  portions  of  the  earth  that  cannot  be  provided 
by  a monostatic  system.  Representative  orbits  and  satellite  con- 
figurations are  analyzed  which  provide  coverage  of  a surveillance 
fence  deployed  along  a constant  latitude.  These  orbits  are  used 
to  generate  bistatic  range-Doppler  maps  of  the  ground  clutter  and 
of  aircraft  on  representative  flight  paths  through  the  surveillance 
fence.  Range-Doppler  maps  of  the  clutter  and  aircraft  form  the 
basis  for  waveform  design.  The  examples  considered  indicate  that 
on  the  order  of  60  dB  of  clutter  rejection  is  needed.  Analog  filter 
techniques  have  been  able  to  achieve  sidelobe  responses  ~40  dB 
below  the  main  lobe  response.  Advances  in  digital  filtering 
techniques  indicate  sidelobe  levels  much  lower  than  40  dB  below 
the  main  lobe  response  of  the  waveform  ambiguity  function  are 
possible.  Examples  of  waveforms  are  presented  which,  with  suitable 
weighting,  could  provide  clutter  rejection  on  the  order  of  the 
needed  60  dB. 
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II. 


SYSTEM  REQUIREMENTS 


The  objectives  of  this  surveillance  system  are  to  search 
for,  detect  and  track  hostile  aircraft  approaching  the  North 
American  continent.  The  primary  approach  corridor  can  be  consid- 
ered to  be  from  the  Soviet  Union  over  the  polar  region  to  North 
America.  Two  additional  attack  corridors  could  also  exist  from 
the  broad  ocean  areas  east  and  west  of  CONUS.  Use  of  these 
corridors  by  the  Soviet  Union  requires  a large  scale  refuel 
operation  over  the  broad  ocean  areas.  Coverage  of  them  is 
considered  a secondary  issue.  They  could  be  covered  at  the 
expense  of  increasing  the  system  capability. 

The  amount  of  warning  time  required  of  a surveillance 
system  is  a function  of  many  variables.  Some  of  these  variables 
are  attacking  aircraft  velocity,  interceptor  velocity,  interceptor 
reaction  time,  minimum  distance  from  CONUS  where  an  intercept 
engagement  may  occur,  etc.  We  chose  arbitrarily  that  the  near 
edge  of  the  surveillance  fence  be  at  least  500  nautical  miles 
from  CONUS. 

2 . 1 Fence  Configuration 

The  surveillance  required  to  meet  the  objectives  of 
searching  for,  detecting  and  tracking  aircraft  approaching  North 
America  could  be  provided  by  many  configurations  of  fences.  This 
analysis  will  attempt  to  evaluate  the  placement  of  a fence  which 
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would  be  capable  of  intercepting  a major  portion  of  threatening 
aircraft  flight  paths  to  North  America.  The  East  - West  (width) 
extent  of  the  fence  depends  on  which  flight  paths  into  North 
America  must  be  detected.  A requirement  also  to  detect  flights 
from  Europe  and  South  America,  for  example,  would  require  a fence 
be  constructed  which  completely  encircles  North  America.  The 
width  of  the  fence  depends  on  aircraft  velocity  and  revisit  search 
time.  The  aistance  between  the  fence  and  the  border  of  North 
America  is  determined  by  the  desired  warning  time  and  aircraft 
velocity. 

The  primary  aircraft  threat  considered  in  this  analysis 
is  an  USSR  long  range  bomber  force.  The  Soviet  airfields  which 
are  suitable  for  and  capable  of  serving  as  launching  and  staging 
areas  for  an  attack  against  North  America  extend  across  the  entire 
continent  of  Asia  from  23®  to  177®  East  Longitude.  A wide  range 
of  attack  azimuths  could  be  employed  against  North  American  targets 
which  places  severe  constraints  on  a fence  configuration  and 
location. 

A bistatic  satellite-borne  radar  system  places  further 
constraints  on  the  type  of  fence  that  can  be  constructed.  The 
satellites  must  be  placed  in  high  altitude  orbits  to  obtain  long 
viewing  times  of  the  fence  with  a minimum  number  of  satellites. 

A low  altitude  orbit  will  require  a much  larger  number  of  satellites, 
which  tends  to  increase  the  complexity  and  cost  but  allows  a 
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graceful  degradation  of  overall  system  performance  if  some 
satellites  are  lost.  In  addition,  a requirement  for  a constant 
bistatic  angle  or  small  range  of  bistatic  angles  as  the  fence  is 
searched  will  place  constraints  on  the  placement  and  configuration 
of  the  surveillance  fence. 

2.2  Typical  Aircraft  Flight  Paths 

The  threat  considered  for  this  analysis  consists  of  USSR 
based  aircraft  attacking  CONUS.  The  threat  may  be  either  a single 
aircraft,  a squadron  or  a massive  coordinated  attack  The  Surveil- 
lance system  should  accommodate  any  size  of  attack.  Without 
attempting  to  establish  likely  paths  between  USSR  airfields  and 
CONUS  strategic  targets,  we  considered  only  great  circle  paths 
from  northern  USf-:  to  CONUS.  Flight  paths  over  the  broad  ocean 
areas  with  in-fl.ght  refueling  are  a possibility.  A system 
designed  to  include  these  cases  will  require  additional  capacity, 
depending  on  the  additional  volume  to  be  searched.  It  is  possible 
that  the  terminal  phase  of  the  flight  might  occur  at  a low  altitude 
and  supersonic  velocities;  however,  this  should  not  affect  our 
analysis  in  any  significant  manner. 

The  aircraft  will  be  assumed  to  fly  at  a nominal  cruise 
altitude,  ( ~10  km)  and  cruise  velocity  (925  km/hr.).  Table  1-1  is 
a matrix  of  flight  paths  chosen  to  represent  typical  attack 
azimuths.  Figure  II-l  shows  the  locations  of  the  targets  and  the 
nominal  origin  of  the  aircraft  (airfield)  with  a straight  line 


TABLE  I-l 


THREAT  SCENARIO 

•> 

'V 


Flight 


Airfield 

Target 

Latitude°N 

Longitude°E 

Latitude°N 

Longitude°E 

55 

157 

35 

-120 

55 

157 

35 

- 75 

65 

75 

40 

-105 

55 

30 

35 

-120 

30 

30 

35 

- 70 

Range 


6475 


8850 


8350 


9610 


7750 


joining  the  target  and  airfield.  Table  1 also  contains  the  great 
circle  distance  between  the  targets  and  airfields. 

2 . 3 Surveillance  Time 

One  of  the  parameters  in  the  radar  equation  which  deter- 
mines system  performance  is  the  solid  angle  search  rate.  A large 
solid  angle  search  rate  requires  large  power-aperture  products. 

In  this  section  we  attempt  to  minimize  the  search  rate  require- 
ments consistent  with  the  objectives  of  the  surveillance  system. 

An  aircraft  flight  path  that  traverses  the  shortest 
distance  through  the  fence  (a  worst  case)  will  determine  the 
maximum  allowable  time  between  successive  looks.  The  aircraft 
is  assumed  to  follow  a great  circle  route  and  does  not  perform 
evasive  maneuvers.  Flight  paths  in  an  easterly  or  westerly 
direction  from  the  USSR  will  most  likely  require  an  in-flight 
refueling,  however,  this  scenario  does  not  impact  the  waveform 
and  clutter  design  requirements.  Additional  radar  fences  could 
be  erected  along  the  satellites  orbital  path  without  increasing 
the  power-aperture  requirements.  Figure  II-2  shows  the  geometry 
of  the  problem  and  defines  nome  terms  for  the  analysis  that  follows. 
The  minimum  angular  extent  of  the  fence  is  0.  The  aircraft  flight 
path  in  the  p' ^ne  cf  Fig.  II-2  represents  a worst  case  condition 
since  all  other  flight  paths  through  the  fence  are  longer,  permit- 
ting longer  search  times. 
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The  altitude  of  the  satellite  is  much  greater  than  the 
altitude  of  the  aircraft.  The  range  from  the  satellite  to  the 
aircraft  can  be  approximated  by  the  range  from  the  satellite  to 
the  earth. 

The  minimum  distance  through  the  fence  is  approximately: 


cos  a 

An  aircraft  traveling  at  a velocity  v^  will  pass  through  the 
fence  in  a time  given  by 

t = — - R9 

V.  ~ V.  cos  a 
t t 

If  we  require  two  independent  looks  at  the  aircraft  while  it  is 
in  the  fence,  the  search  period  becomes 

R9 

^s  * 2v^  cos  a 

The  maximum  allowable  search  time  to  insure  at  least  two  looks 
at  the  aircraft  is  plotted  in  Fig.  II-3  vs  the  distance  the  air- 
craft must  travel  through  the  fence.  The  velocity  of  the  aircraft 
is  a parameter.  As  the  distance  through  the  fence  increases,  the 
allowable  search  time  increases  proportionally.  If,  for  example. 
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3000 


Fig.  II-3. 


FENCE  DISTANCE  ^ (km) 
Searcli  time  requirements. 
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a 925  km  (500  nautical  mile)  fence  were  used  and  the  average 
aircraft  velocity  were  925  km/hour,  the  fence  would  have  to  be 
searched  every  half  hour. 

2 . 4 Potential  Advantages  of  Bistatic  Systems 

The  principal  potential  advantages  of  a bistatic  radar 
for  aircraft  detection  are  aircraft  cross-section  enhancement 
and  ground  clutter  cross-section  density  diminution.  The  combina- 
tion of  these  two  potential  effects  is  to  increase  the  signal-to- 
clutter  ratio,  S/C,  and  thereby  simplify  the  detection  process. 

The  existence  and  extent  of  these  potential  advantages  must 
ultimately  be  determined  by  measurements. 

Bistatic  radar  cross-soctions  of  simple  shapes  have  been 
calculated  in  the  literature. The  forward  scatter  cross-section 
of  simple  shapes  has  been  shown  to  be  many  times  the  backscatter 
cross-section  when  the  wavelength  is  small  compared  to  the  target 
dimensions.  However,  it  is  difficult  to  predict  whether  a compli- 
cated target,  such  as  an  aircraft,  will  exhibit  this  large  cross- 
section  enhancement. 

(2) 

Bistatic  clutter  measurements  have  been  made  for 
various  terrain  conditions  but  only  over  limited  bistatic  geometries. 
Additional  data  are  needed  before  definitive  conclusions  about 
potential  advantages  of  bistatic  clutter  cross-section  can  be  made. 


12 


III.  RADAR  PARAMETER  ESTIMATES 

In  this  section,  the  size  of  a radar  to  perform  search 
from  a satellite-borne  platform  is  estimated  as  a function  of 
search  rate.  Several  example  designs  are  presented  and  some 
comments  about  system  realizability  are  made.  Included  are  some 
remarks  about  communications  and  control.  The  results  of  the 
power-aperture  product  estimates  are  parameterized  over  a large 
range  intej-val  to  include  satellite  orbits  up  to  geosynchronous 
altitude.  At  geosynchronous  altitude,  the  ground  clutter  will 
appear  at  zero  Doppler  with  a small  spread.  In  general,  the 
target  will  present  a Doppler  different  from  the  clutter  Doppler. 
At  Ijwer  satellite  altitudes,  the  clutter  Doppler  will  be  offset 
from  zero  frequency  and  will  be  a function  of  the  angle  the  radar 
line  of  sight  makes  with  the  velocity  vector.  This  tends  to 
complicate  Doppler  filtering.  Furthermore,  a geostationary 
satellite  radar  system  allows  a fence  to  be  placed  at  any  desired 
location  within  the  field  of  view  and  its  location  can  be  changed 
at  any  time.  It  provides  for  continuous  viewing  of  a specific 
area  and  allows  for  continuous  tracking  of  detected  targets  for 
extended  periods  of  time.  In  addition  it  allows  for  repeated 
looking  at  specific  areas  such  as  most  likely  flight  paths.  The 
main  disadvantage  of  a geosynchronous  altitude  system  is  that  a 
very  large  power-aperture  product  is  required,  as  will  be  seen 
in  this  section.  A lower  altitude  system  requires  a lower  power- 
aperture  product;  however,  more  satellites  are  requiiod  to  provide 
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the  desired  coverage. 

3.1  Search  Power-Aperture  Product 

A useful  form  of  the  radar  equation  for  sizing  a search 
radar  system  operating  in  a bistatic  mode  is  given  by 


PA 


4tt  k T L 


(1) 


where 

n a:  solid  angle  search  rate 

Rj.  = range  from  target  to  receiver 

= range  from  target  to  transmitter 

s 

^ =«  signal  to  noise  ratio 

L = system  losses 

T ® system  temperature 

k = Boltzmann's  constant 

= bistatic  target  cross-section 

A = aperture  area 

P = average  radiated  power 

The  equation  can  be  put  in  a more  useful  form  if  we  relate  the 
solid  angle  search  rate  to  the  ground  coverage  rate.  The  inter- 
section of  the  beam  with  the  surface  of  the  earth  determines  the 


ground  coverage  per  beam  position.  For  a narrow  circular  beam, 
this  area  is  given  by 

2 2 

Area  w (2) 

cos  a 

where 

a = angle  the  radar  LOS  makes  with  the  target 
local  vertical 

2 

9 = solid  angle  subtended  by  tne  beam 

(omitting  a factor  »it/4) 

R = range  to  target 

For  a bistatic  system  the  area  illuminated  is  determined  by  the 
smaller  of  the  transmitter  or  receiver  illumination.  If  the 
transmitter  and  receiver  are  at  the  same  altitude,  hence  at  nearly 
equal  ranges  in  most  cases,  the  smaller  beamwidth  will  determine 
the  illuminated  area.  Differentiating  both  sides  of  equation  (2) 
and  assuming  the  range  and  radar  angle,  a , to  remain  essentially 
constant  over  the  search  volume,  we  obtain 

d Area  ^ ^ d (9^) 

dt  cos  a * dt 

Substituting  equation  (3)  into  equation  (1)  we  have 


15 


_ 4Tr  k TL  ,S^  „2  /d  Area  \ cos  a (4) 

« \-^j 


An  aircraft  attack  against  CONUS  would  most  likely  use 
routes  across  the  northern  portion  of  Canada.  A small  portion  of 
these  routes  might  occur  over  the  broad  ocean  areas;  however,  for 
this  analysis  only  flights  over  Northern  Canada  are  considered. 

The  angle  the  radar  line  of  sight  (from  geostationary  altitude) 
makes  with  the  target  local  vertical  '^aries  between  64  and  82 
degrees  for  a 100  n.mi  wide  fence  centered  about  65°  N latitude. 
Therefore,  an  average  value  for  a of  70  degrees  will  be  used.  The 
ground  coverage  search  rate  ( — — ) is  plotted  in  Fig.  III-l 
as  a function  of  the  power  aperture  product  given  by  equation  (4) 
with  range  a parameter. 

The  strong  dependence  of  the  power-aperture  product  on 
search  rate  is  demonstrated  in  Fig.  III-l  and  illustrates  the  need 
to  minimize  the  search  rate  to  maintain  system  realizability.  The 
ground  search  rate  is  given  by; 


d Area 
dt 


T X W 
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te  vs.  power  aperture  product. 
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where 


T 

W 

t. 


thickness  of  fence 

width  of  fence 

time  to  search  entire  fence 


T 


n V, 


where 


n 

'^t 

d Area 
dt 


number  of  opportunities  provided  to  observe 
a target  (uniform  search  assumed) 

target  velocity 


n Vt  W 


Note  that  the  ground  search  rate  is  independent  of  fence  thickness. 
It  depends  only  on  the  number  of  independent  opportunities  to  view 
the  target  (assuming  a uniform  search) , the  aircraft  velocity  and 
the  width  of  fence.  For  this  analysis,  we  will  assume  an  average 
aircraft  velocity  of  500  n.  mi. /hour  (0.26  km/s) . If  we  require 
at  least  two  opportunities  to  view  the  aircraft  while  it  is  in  the 
fence  and  we  assume  that  the  fence  will  be  covered  once  per  hour, 
the  fence  must  be  1000  n.  mi.  (1800  km)  thick.  The  width  of  the 
fence  should  be  such  that  it  intercepts  the  major  portion  of 
possible  direct  flight  paths  fiom  the  Soviet  Union  to  the  United 
States.  This  implies  a fence  approximately  3000  n.  mi.  (5509  km) 
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wide.  The  ground  search  rate  for  this  set  of  conditions  is 

6 3 

3 X 10  sq.  n.  mi. /hour  or  2.85  x 10  sq.  km/s.  Note  that 

wavelength  does  not  appear  explicitly  in  equation  (4)  so  that 

this  result  applies  at  any  wavelength  provided  the  system 

temperature,  losses  and  target  cross-section  used  are  the  same  at 

these  wavelengths.  The  wavelength  dependence  is  implicit  in  the 

expression  for  aperture. 
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The  wavelength  dependence  of  the  power  and/or  aperture 
can  be  plotted  once  an  operating  altitude  (range)  or  a beamwidth 
is  chosen.  Figures  III-2  and  III-3  are  example  plots  showing  the 
power  and  aperture  dependence  on  wavelength  at  nominal  ranges  of 
5000  and  40000  km  respectively.  Antenna  beamwidth  is  a parameter. 
The  power-aperture  product  increases  by  18  dB  when  the  range  is 
increased  from  5000  km  to  geosynchronous  altitude  ranges. 

3. 2 Track  Power-Aperture  Product 

The  requirements  for  a tracking  radar  system  are  described 
by  the  following  form  of  the  radar  equation 


P^A- 


4tt  k T L B „ 2 „ 2 ,2 


(6) 


I 
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For 


and 


B = 


1 

T 


= Pt 


let  PRF  = BRF  (Burst  repetition  frequency)  equation  (6)  becomes 

PA^  = 4jl.|  T JL  ^2 

b 

A pulse  burst  waveform  is  needed  for  clutter  filtering  and  is 
discussed  in  Section  VI.  The  trade  rate  is  1 burst  per  dwell. 

If  the  bistatic  transmitter  and  receiver  are  in  satellites  at  the 
same  altitude,  the  ranges  from  the  transmitter  to  :he  target  and 
the  receiver  to  the  target  will  be  nearly  equal  in  most  cases. 

For  this  analysis,  we  set  = Rj.  = R.  The  aperture  and  wavelength 
are  related  to  the  beamwidth  by  equation  (5) . Substituting  this 
relationship  into  ecuation  (7)  and  rearranging  terms  we  have 


PA 


4lJ^  (|)  r' 


(BRF) 


(8) 


Figure  III-4  plots  the  burst  repetition  frequency,  BRF,  as  a 

2 

function  of  the  tracking  performance  parameter  PA/9  with  range 
a parameter. 
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The  track  function  requirement  is  in  addition  to  the 
primary  surveillance  requirement.  The  additional  system  capability 
to  perform  the  track  function  can  be  easily  determined  using 
Fig.  III-4  and  the  search  radar  beamwidth.  The  revisit  time,  or 
1/BRF,  is  determined  by  the  target  velocx^y,  antenna  beamwidth 
and  range  to  the  target.  The  most  stringent  requirement  is  when 
the  target  is  allowed  to  maneuver  between  pulse  bursts.  In  this 
case,  the  BRF  becomes; 


If  the  target  is  restricted  to  fly  a great  circle  route,  the  BRF 
track  requirement  becomes  a function  of  the  tracking  algorithm  and 
the  measurement  accuracy.  In  general,  this  requirement  will  not 
be  as  severe  as  that  given  for  the  maneuvering  target  example 
above . 

As  an  example,  consider  the  radar  system  at  geosynchronous 

altitude.  The  target  is  assumed  to  be  moving  at  500  n.  mi. /hour 

-3  -1 

and  the  beamwidth  is  0.1  degree.  The  BRF  becomes  3.68  x 10  sec  . 
If  we  allow  the  target  to  move  a maxinium  of  one-half  beamwidth 
between  successive  track  pulses,  the  average  power  becomes  30.4  watts 
per  target.  An  operating  wavelength  of  0.1  meter  was  assumed. 

This  represents  a small  increase  in  average  power;  however,  if 
many  targets  (>100)  must  be  tracked,  a sizeable  increase  in  system 
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performance  will  be  necessary. 


3.3  System  Realizability 

The  power-aperture  curves  presented  in  the  previous 
sections  form  the  basis  for  determining  the  tradeoff  between  power, 
aperture  and  range.  The  other  parameters  in  the  search  radar 
equation  are  fairly  well  defined  constants.  The  parameter  least 
defined  is  the  bistatic  radar  cross-section.  It  remains  to  be 
better  defined  by  measurements.  In  this  analysis  we  have  assumed 
the  target  cross-section  is  15  dBsm  and  remains  constant  with 
changes  in  aspect  angle  and  wavelength. 

A satellite-borne  radar  system  viewing  targets  against  a 

background  of  the  earth  will  have  to  detect  targets  imbedded  in 

clutter.  The  amount  of  clutter  competing  with  the  target  is 

determined  by  the  resolution  cell  size  in  range  and  cross-range 

(Doppler)  and  the  reflectivity  (o^)  of  the  clutter.  Measurements 
(3) 

have  been  made  of  isolated  samples  of  ground  clutter  at  L-band 
and  X-band.  These  data  are  for  small  bistatic  angles  and  are 
mainly  in  the  plane  containing  the  transmitter  and  local  vertical. 
Wide  bistatic  angle  data  out  of  the  plane  described  above  does 
not  exist  so  far  as  we  have  been  able  to  determine.  The  existing 
data  suggest  that  there  is  a large  variation  in  the  bistatic  clutter 
cross-section  and  that  there  is  little  systematic  difference  between 
it  and  monostatic  clutter  cross-section.  A satellite-borne  bistatic 
system  will  be  viewing  large  areas  which  are  composites  of  various 
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types  of  surfaces.  It  is  difficult  to  predict  the  behaviour  of 
the  bistatic  clutter  return  as  the  satellite  passes  over  varying 
terrain  such  as  cities,  smooth  fields,  mountains,  ice,  snow,  etc. 

The  available  data  suggest  that  the  resolution  cell  size 
should  be  made  as  small  as  practicable.  The  techniques  for  achiev- 
ing resolution  in  the  range  direction  are  well  known.  Bistatic 
range  resolution  on  the  order  of  the  maximum  dimensions  of  the 
target  (30m  < 6 < 100m)  is  well  within  the  state  of  the  art.  The 
resolution  in  the  cross-range  direction  is  obtained  by  making  the 
antenna  beamwidth  small.  This  implies  a large  aperture  which  may 
not  be  easily  achieved  on  a spaceborne  platform.  The  largest 
antenna  known  to  us  to  be  deployed  aboard  a space  platform  is  a 
10  meter  diameter  parabolic  dish  structure.  Techniques  which  are 
tolerant  of  mechanical  errors  and  still  capable  of  maintaining  a 

reasonable  antenna  beamwidth  (gain)  will  have  to  be  exploited  for 

(4) 

space-borne  applications. 

The  intersection  of  the  antenna  beam  and  the  surface  of 
the  earth  will  be  approximately  an  ellipse  with  the  major  axis 
oriented  along  the  range  direction  and  the  minor  axis  normal  to 
the  range  direction.  The  cross-range  dimension  vs  range  to  the 
satellite  is  plotted  in  Fig.  III-5  with  antenna  beamwidth  a 
parameter.  At  geosynchronous  altitudes  (range  = 40,000  km)  , an 
antenna  beamwidth  of  0.1  degree  gives  a cross-range  dimension  of 
approximately  70  km. 
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Fig.  III-5. 


Cross  range  vs.  range  to  satellite. 


The  parameters  used  to  generate  the  power-aperture  vs 
wavelength  curves  of  the  previous  section  will  be  used  in  a sample 
design  for  a bistatic  radar  system  at  geosynchronous  altitude. 

The  power-aperture  product  to  search  a 1000  x 3000  n.  mi.  ground 
area  once  per  hour  is 


PA  = 2.7  X 10^  w m^ 

For  this  example,  let  us  assume  the  following  additional  parameters. 

2 

= clutter  reflectivity  = -10  dBsm/m 

e = beamwidth  = 0,1  degree 

R = range  to  target  = 40,000  km 

6j,  = range  resolution  = 100  m 

s 

^ = required  target-to-clutter  ratio  after 

filtering  = +15  dB 

The  bistatio  clutter  cross-section  in  a range  cell  can  be  approxi- 
mated, for  small  bistatic  angles,  by 

0 ~ a A = o (RO)  (6  r) 

CO  o 

o ~ 60  dBsm 
c 
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To  achieve  a target  to  clutter  ratio  of  +15  dB  on  a target  cross 
section  of  15  dBsm,  a clutter  filter  with  a -60  dB  sidelobe 
response  is  required.  Other  examples  are  discussed  in  Section  V. 

Let  us  arbitrarily  choose  an  operating  frequency  of  S-band 
(wavelength  = 0.1  m) . The  aperture  to  generate  a 0.1  degree  circular 
beamwidth  becomes 

2 

A = ^ = 3282  m^ 

0 

The  required  average  power  is 
P » 825  watts 

The  system  parameters  are  summarized  in  Table  III-l.  The  para- 
meters of  this  sample  system  are  ambitious.  The  generation  of 
approximately  1 kw  average  power  for  an  extended  period  of  time 
in  space,  deployment  of  a large  aperture  in  space  and  advanced 
clutter  filtering  techniques  would  need  to  be  developed.  Some 
filter  designs  are  discussed  in  Section  VI. 


TABLE  II 1-1 


An  Example 

Radar  System  Parameters 


Average  Power 

Wavelength 

Beam  Width 

Antenna  Diameter 

Range  Resolution 

System  Temperature 

Losses  (Transmit  & Receive) 

Target  cross-section 

Range 

Search  Time 
Search  Area 

Clutter  Filter  Response 

S/C  required  after  filtering 

S/N 


825  watts 
0.1  meter 
0.1  degree 
60  m 
100  m 
1000°  K 
10  dB 
15  dBm^ 

40,000  km 
1 hour 

1.03  X lO’^'^  m‘ 

;^60  dB  clutter  rejection 
+ 15  dB 
+ 15  dB 


IV. 


REPRESENTATIVE  SATELLITE  ORBITS 


4.1  General  Categories 

Satellite  orbits  for  bistatic  surveillance  can  be  grouped 
into  several  categories  which  relate  to  various  system  design 
objectives.  Low  altitude  satellites  allow  shorter  ranges  and 
therefore  reduced  power-aperture  product  but  at  the  expense  of 
small  coverage  per  satellite.  Geostationary  satellites  are  always 
in  position  to  scan  a fence  but  their  high  altitude  results  in 
large  power-aperture  requirements.  Intermediate  altitude  satellites 
may  provide  a compromise  between  coverage  and  power-aperture 
requirements.  Satellites  whose  periods  are  commensurate  with  that 
of  the  earth  will  have  repetitive  ground  traces  and  coverage.  This 
characteristic  is  desirable  to  minimize  the  number  of  satellites, 
but  not  essential.  Additional  objectives  are  to  select  bistatic 
geometries  which:  (1)  ease  the  clutter  rejection  problem  by 

advantageous  placement  of  the  clutter  relative  to  targets  in  the 
range-Doppler  (R-R)  plane  and  (2)  exploit  any  potential  target 
bistatic  cross-section  enhancement  and/or  clutter  diminution. 

4. 2 Virtual  Radar  Location 

The  location  of  a satellite  borne  bistatic  radar  can  be 
characterized  by  the  sub-satellite  point  of  a virtual  satellite- 
borne  monostatic  radar  on  the  bisector  of  the  lines-of-sight  to 
an  arbitrarily  located  target.  The  range  of  the  virtual  radar 

% 
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from  the  target  is  selected  so  that,  for  all  other  parameters 
of  the  bistatic  and  virtual  monostatic  radars  the  same,  the  S/N 
ratio  is  the  same  for  both  radars.  This  range  is  the  geometric 
mean  of  the  receiver-target  and  transmitter-target  ranges  of  the 
bistatic  radar.  The  virtual  radar  sub-satellite  point  is  illustrated 
in  Fig.  IV-1. 

A second  characterization  of  the  virtual  radar  location 
is  the  sequence  of  points  on  the  earth's  surface  at  which  the 
clutter  can  be  separated  from  an  aircraft  by  range  gating.  The 
local  vertical  at  these  points  lies  along  the  bistatic  range 
vector  and  their  locus  is  generated  by  the  motion  of  the  satellites 
and  the  earth.  Figure  IV-2,  illustrates  the  geometry  for  one  of 
these  points.  For  a satellite-borne  monostatic  radar  these  points 
are  the  sub-satellite  points. 

Examples  of  coverage  from  bistatic  radars  are  given  below 
for  representative  orbits.  The  R-R  characteristics  of  clutter  and 
aircraft  are  discussed  in  Section  V. 

4.3  Examples  of  Coverage 

The  virtual  radar  location  as  defined  above  provides  a 
convenient  means  of  indicating  the  viewing  aspects  of  various 
points  on  the  earth's  surface.  Plots  of  these  locations  for  four 
representative  bistatic  configurations  are  presented  in  Figs.  IV-3 
to  IV-9. 
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18-3-19098 


VIRTUAL 

RADAR 


Fig.  IV-2,  Geometry  illustrating  coincidence  of  virtual  radar 
sub-point  and  observed  point. 


Figure  IV-3  shows  the  sequence  of  virtual  radar  locations 
of  a bistatic  radar  viewing  five  points  at  60°  N latitude  and 
within  the  surveillance  fence.  The  receiver  and  transmitter  of 
the  bistatic  radar  are  each  on  a 12  hour  period  satellite  whose 
elements  are  included  in  Fig.  IV-3.  Note  that  the  nodes  are  180° 
apart.  The  periods  and  inclinations  of  these  j .tellites  are  equal 
and  are  adjusted  so  that  their  angular  rates  at  apogee  match  the 
angular  rate  of  the  Earth.  The  periods  are  also  made  commensurate 
with  that  of  the  Earth.  Because  of  this  matching  of  angular  rates 
the  virtual  radar  location  does  not  move  greatly  during  the  time 
the  viewed  points  are  visible  to  the  satellites.  Also,  since  the 
satellites  are  on  opposite  sides  of  the  North  Pole,  the  bistatic 
range  vector  to  the  pole  lies  along  the  Earth's  axis  and  consequently 
the  pole  is  visible  (about  9 hours) , that  is  the  virtual  radar  is 
stationary  over  the  pole. 

r.qure  IV-4  similarly  shows  the  virtual  radar  locations  for 
two  6 hr.  period  satellites  in  the  same  plane,  nodes  180°  apart, 
moving  in  opposite  directions  and  with  apogees  at  105°  E longitude. 
The  virtual  locations  tend  to  be  at  a constant  latitude  (the  orbit 
plane  inclination)  during  much  of  the  time  the  viewed  points  are 
visible.  Figure  IV- 5,  for  the  same  pair  of  satellites,  shows  the 
sequence  of  locations  for  which  clutter  can  be  range-gated.  This 
sequence  is  nearly  at  constant  latitude  and  covers  about  half  of 
•.  le  fence  width  in  4.4  hours.  This  is  a consequence  of  the  virtual 
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18-3-19099 


Fig.  IV-3.  Ground  traces  of  virtual  satellite  observing  various 
points:  example  1. 


LATITUDE  (deg) 
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Fig.  IV-4.  Ground  traces  of  virtual  satellite  observing  various 
points:  example  2. 
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BISTATIC  SATELLITES 
1 2 


APOGEE  HEIGHT  (km) 

19,906 

19,906 

PERIGEE  HEIGHT (km) 

400 

400 

INCLINATION  (deg) 

65 

115 

ASCENDING  NODE  (deg) 

- 1 57.4 

22. 

ARGUMENT  PERIGEE  (deg) 

270 

270 

PERIOD  (hr) 

6 

6 

TIME  TICKS  *0.1  hr 
TOTAL  TIME  • 4.4hr 


LONGITUDE  (deg) 

Fig.  IV-5.  Sequence  of  points  for  which  Clutter  can  be  Range 
Gated. 
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radar  being  nearly  fixed  in  inertial  space  for  this  example. 

Although  the  location  of  this  sequence  is  within  the  fence,  many 
satellites  would  be  required  to  provide  similar  virtual  radar  traces 
to  fill  out  the  complete  width  and  length  of  the  fence.  It  is 
unlikely  that  the  virtual  radar  location  of  either  Fig.  IV-4  or 
Fig.  IV-5  can  be  closely  approximated  by  a monostatic  radar, 
although  operationally  equivalent  performance  may  be  possible  by 
altering  the  shape  of  the  fence.  The  bistatic  angle  for  this 
system  will  vary  with  time  from  0°  to  near  180°. 

Figures  IV-6  and  IV-7  show  virtual  radar  locations  for  a 
geostationary  satellite  at  -105°  E longitude  and  a 6 hr.  period 
polar  satellite  that  has  its  apogee  over  the  pole.  The  traces  of 
Fig.  IV-6  essentially  overlap  so  that  the  virtual  radars  for  all 
of  the  five  points  in  the  fence  have  comparable  "orbits”.  The 
trace  of  Fig.  IV-7  indicating  points  where  clutter  may  be  range 
gated  is  mostly  at  low  latitudes  and  hence  not  useful  as  a fence. 

Figures  IV-8  and  IV-9  are  for  a pair  of  6 hr.  period 
satellites  with  65°  inclination  but  with  nodes  90°  apart.  These 
traces  are  similar  to  those  of  Figs.  IV-6  and  IV-7  but  have  greater 
extent  in  both  latitude  and  longitude.  The  trace  of  Fig.  IV-9  is 
suggestive  of  a fence  covering  East  and  West  approaches  to  the  US 
as  well  as  from  the  North. 

The  virtual  radar  locations  for  low-altitude  co-orbital 
satellites  lie  on  a narrow  strip  above  the  ground  trace  of  the 
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-180  -160  -140  -120  -100  -80  -60  - 40 

LONGITUDE  (deg) 

Fig.  IV-6.  Ground  traces  of  a fictitious  satellites  with  a 
monostatic  radar  equivalent  to  a bistatic  radar  observing 
various  points. 


TUDE  (deg) 


satellites  and  are  not  very  different  from  the  orbit  of  a monostatic 
radar  midway  between  the  two  satellites  (and  somewhat  higher) . As 
the  spacing  of  the  co-orbital  satellites  is  reduced,  the  virtual 
radar  location  converges  to  the  midpoint  between  the  two. 

These  examples  indicate  that  virtual  radar  locations  can 
be  established  at  various  positions  relative  to  viewed  points  in 
the  fence  and  that  a high  altitude  virtual  radar  can  be  made  to 
remain  "on  station"  in  rather  narrow  regions  during  much  of  the 
time  the  points  are  visible.  Lines  along  which  clutter  can  be 
range-gated  can  also  be  favorably  located.  The  possibility  of 
broadening  these  lines  into  a swath  for  targets  at  some  altitude 
above  the  Earth,  and  for  viewing  large  segments  of  the  line 
instantaneously  should  be  investigated.  A sufficiently  wide  swath 
and  large  instantaneous  line  segment  coverage  may  permit  a reason- 
able number  of  satellite  radars  to  scan  a fence  using  range-gating 
only  to  eliminate  the  clutter.  (Appendix  B addresses  this  technique 
for  monostatic  radars.) 

The  requirements  for  beam  steering  also  influence  orbit 
selection.  Since  precise  scanning  of  large  antennas  in  space 
does  not  eppear  attractive  it  is  expected  that  the  antenna  would 
be  stabilized  and  electronically  scanned.  To  keep  the  effective 
aperture  large  the  scan  angle  from  boresight  should  be  kept  small 
(a  60°  scan  angle  from  boresight  approximately  results  in  a 3 dB 
aperture  loss) . The  fence  scan  angle  and  corresponding  aperture 
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losses  can  easily  be  made  sufficiently  small  (less  than  0.6  dB) 
by  making  the  satellite  altitudes  at  least  twice  the  longest 
dimension  of  the  fence.  However,  as  the  satellites  move  toward 
the  horizon  of  the  fence,  the  of f-boresight  angles  increase. 
Consequently,  to  maintain  low  aperture  loss  the  satellite  altitudes 
must  be  increased  or  the  useful  viewing  time  of  the  fence  reduced. 
An  inertially  stabilized  antenna  in  a geostationary  orbit  appears 
best  suited  for  simple  and  efficient  aperture  use. 


45 


V. 


RANGE-DOPPLER  MAPS  AND  WAVEFORM  DESIGNS 


The  satellite-borne  bistatic  aircraft  surveillance  radar 
searches  for  aircraft  against  an  earth  background  and  therefore 
must  be  able  to  distinguish  aircraft  from  background  clutter.  The 
aircraft  signal  may  be  within  the  range  and  Doppler  extent  of  the 
clutter  or  may  be  separate  from  it  in  either  range  or  Doppler  or 
both.  Range-Doppler  or  R-R  maps  of  the  clutter  and  aircraft  signals 
can  be  made,  to  show  their  relationship.  The  maps  can  then  be  used 
to  select  suitable  clutter  rejection  techniques.  A description  of 
the  computation  of  bistatic  range  and  Doppler  for  making  R-R  maps 
is  given  in  Appendix  A. 

There  is  the  possibility  of  range-gating  the  clutter  for 
certain  viewing  geometries  and  aircraft  altitudes.  A special  case 
of  this  was  discussed  in  Section  IV  for  the  virtual  monostatic 
radar  viewing  the  earth's  surface  at  normal  incidence.  A discussion 
of  clutter  range  gating  at  oblique  viewing  angles  for  a monostatic 
radar  is  given  in  Appendix  B and  illustrates  concepts  that  can  be 
extended  to  bistatic  radars. 

The  primary  technique  considered  here  for  separation  of 
aircraft  and  clutter  will  be  some  form  of  Doppler  filtering.  The 
effectiveness  of  waveforms  will  depend  on  the  distribution  of  clutter 
in  Doppler,  the  location  of  aircraft  in  Doppler  and  signal  integra- 
tion time.  The  clutter  Doppler  characteristics  are  dependent  on 


46 


I 


viewing  geometry  and  beam  parameters  whereas  the  aircraft  Doppler 
is  only  geometry  dependent. 


The  waveforms  sought  must  not  only  provide  the  needed 
clutter  reduction  but  must  also  be  consistent  with  the  required 
search  rate  for  the  fence  as  derived  in  Section  III.  That  is  the 
illuminated  area  and  integration  or  dwell  time  must  satisfy  the 
coverage  rate  requirement 


dA 

dt 


Illuminated  area  ^ , 2, 

- Dwai  Hii i 2.9x10  km  /s 


if  a single  radar  is  to  be  able  to  conduct  the  search. 


Figure  V-1  illustrates  the  geometry  for  bistatic  radars 
with  both  fixed  and  time  varying  Doppler  spread  of  the  clutter. 

The  clutter  Doppler  distribution  of  a bistatic  radar  on  low-altitude 
circular  orbits  is  fixed,  dependent  only  on  beam  dimensions.  As 
different  terrain  passes  through  the  f ield-of-view  there  will  be 
large  changes  in  clutter  cross-section,  however,  changes  in  the 
clutter  spectrum  are  expected  to  be  rimall  compared  to  the  Doppler 
spread  due  to  satellite  motion  relative  to  the  earth.  The  clutter 
spectrum  of  a geostationary  bistatic  radar  is  fixed,  centered  at 
zero  Doppler  and  has  zero  Doppler  spread.  Varying  Doppler  spreads 
occur  when  the  beams  are  scanned  or  when  the  satellites  are  on 
elliptic  orbits,  and  are  latitude  dependent  when  the  satellites' 
rates  are  comparable  to  that  of  the  earth.  Varying  clutter 
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distributions  may  require  a time  varying  filter  to  adequately 
reject  the  clutter. 

5.1  Doppler  Resolution  Limitations 

Three  characteristic  relations  exist  between  the  aircraft 
and  clutter  Doppler  i.i  a given  range  cell;  1) aircraft  Doppler 
widely  separated  from  the  clutter;  2) aircraft  Doppler  separated 
but  close  to  that  of  the  clutter;  3) aircraft  Doppler  within  the 
clutter  Doppler  distribution.  When  the  clutter  and  aircraft  are 
separated  in  Doppler  the  Doppler  resolution  of  the  waveform  and 
filter  need  not  be  great  since  a highly  precise  determination  of 
the  velocity  of  the  target  is  not  necessary.  However,  low 
sidelobes  are  essential.  If  the  target  and  clutter  are  closely 
spaced  in  Doppler  the  waveform  and  clutter  filter  must  also  have 
sharp  skirts.  The  coherent  integration  time  will  limit  both  the 
Doppler  resolution  and  filter  skirt  sharpness. 

When  the  aircraft  and  clutter  Doppler  overlap  then  Doppler 
resolution  can  increase  the  S/C,  but  only  by  the  ratio  of  the  total 
clutter  Doppler  spread  to  the  Doppler  spread  of  the  aircraft.  For 

a m^nosatic  radar  viewing  a bomber  with  a characteristic  dimension 

( 5 i 

of  30  m,  the  RMS  Doppler  error  due  to  scintillation  is  estimated 
to  be 


2.72  ( 


2iie . 
A ’ 


n 

'360 


0.56/A  Hz 
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where  e is  the  RMS  angle  noise  (in  linear  units) ~ 0.3  L = 9 m 

fi  is  the  RMS  angular  velocity  ~ 1.3°/s 
A is  the  wavelength,  m 
L is  the  bomber  length  ~ 30  m 

In  terms  of  Doppler  velocity 

AR  'v  0. 3 m/s. 
s 

For  a bistatic  system  it  can  be  expected  that 
ARsb  0 . 3 cos  3/2  m/s 

where  3 is  the  bistatic  angle.  The  clutter  Doppler  spread  is  on 
the  order  of 

ARc  Vgfl  cos  3/2 

where  v^  is  the  velocity  of  the  satellites  relative  tv.  the  earth 
in  m/s  and  0 is  the  beamwidth  in  radians. 

The  S/c  improvement  by  Doppler  filtering  for  this  case  is  then 
given  by 


AR 


c 

sb 


% 


0.  3 
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In  order  to  make  this  ratio  large  one  wants  to  make  large 

by  using  low-altitude  satellites  and  0 large.  However,  to  minimize 

the  total  cross-section  6 should  be  as  small  as  possible.  For 

example,  if  one  used  a 0.6^  beam  (0.01  radian)  and  a 200  km 

altitude  circular  orbit  satellite  (v  = 7700  m/s)  then  the  Doppler 

s 

filtering  improvement  becomes 

V ^ iI?00H0-01)  = 257  24  dB 

0.3  0.3 

A typical  improvement  needed  is  60  dB.  Consequently,  it  appears 
that  aircraft  detection  is  possible  only  if  the  aircraft  Doppler 
does  not  overlap  that  of  the-  clutter. 

5.2  Range  Resolution  Limitations 

The  clutter  competing  with  the  aircraft  is  the  total 
clutter  within  the  illuminating  beam  reduced  by  the  ratio  of  the 
bistatic  range  resolution  to  the  bistatic  range  extent  of  the  clutter. 
However,  the  range  cell  cannot  be  made  arbitrarily  small.  If  it  is 
smaller  than  the  aircraft  dimensions  then  the  aircraft  cross-section 
will  decrease  in  approximately  the  same  ratio  as  the  clutter  cross- 
section  and  there  will  be  no  signal/clutter  improvement.  If  there 
is  signal  integration,  then  the  range  cell  must  encompass  not  only 
the  aircraft  but  the  distance  it  travels  during  the  integration  tune. 
If  the  integration  time  is  less  than  about  0.01  s (to  assure  an 
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adequately  high  ground  coverage  search  rate)  then  for  a typical 
30  meter  length  aircraft  and  300  m/s  velocity  the  range  cell  must 
encompass  a distance  Ax  ^ 30  cos  <p , where  <p  is  the  angle  between 
the  aircraft  velocity  and  the  bistatic  direction.  (Because  of 
aircraft  thickness  Ax  ^ 5 m.)  The  bistatic  range  resolution 
corresponding  to  a length  x along  the  bistatic  range  direction  is 

ARg  + Ax  cos  8/2  ~ 30  cos  ((>  cos  8/2 

where  8 is  the  bistatic  angle.  At  large  bistatic  angles  AR_ 

O 

required  for  a given  Ax  can  become  very  small  for  a given  Ax  so 
that  large  bandwidth  waveforms  would  be  required.  For  maximum 
clutter  reduction  the  range  resolution  could  be  made  to  vary  with 
changes  in  4>  and  8.  However,  choosing  ARg  30  m is  a good  estimate 
for  modest  values  of  (t>  and  8. 

5. 3 Bistatic  Clutter  Cross~Section 

The  bistatic  clutter  cross-section,  a , for  an  elemental 

c 

ground  area,  dA,  is  given  by 

o„  - 0 dA 

c 


or 


0 

c 


Y cos  4>  dA 
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where  a°  is  the  bistatic  cross-section  density  of  the  illuminated 
area,  y is  the  cross-section  density  of  the  area  projected  on  the 
bistatic  range  and  cf)  the  bistatic  angle  of  incidence.  Both  a°  and 
Y are  functions  of  the  geometry,  polarization  and  terrain. 

The  bistatic  clutter  density  in  the  R-R  plane  can  be 
determined  for  specific  bistatic  geometries  from  the  Jacobian  of 
the  transformation  from  ground  coordinates  to  R-R  coordinates,  if 
J is  the  Jacobian  then  the  clutter  cross-section  is 

a = 0°  dA  = 0°  d X dy  = 0°  J dR  dR 

c 


or 


0 = Y cos  ({!  dA  = Y cos  fj)  J dR  dR 

Since  both  J and  cos  41  are  only  geometry  dependent  they  can  be 
combined  so  that  the  latter  form  becomes 


o - y J,  dR  dR 
c ' <p 


J and  J can  be  considered  relative  cross-section  densities  since 


0°  dR  dR 

*See  Appendix  A for  a derivation  of  the  Jacobian. 
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% f « • > 


« •* 


J 


and 


Y dR  dR 

The  unit  of  J and  J,  is  seconds. 

<P 

Typically  is  less  than  -10  dB^^^  for  large  bistatic 
angles  of  incidence  (low  grazing  angles)  and  y ^ o°  since  cos  c|)  ;<  1. 
There  may  be  a strong  specular  when  the  normal  to  the  clutter  is 
along  the  bistatic  range  direction.  Consequently,  the  clutter 
from  a patch  of  ground  with  dimensions  comparable  to  that  of  an 
aircraft  is 


= a°  dA  -V  lo"^  (30  x 30)  « 90  ^ 20  dBsm 

Clearly  an  aircraft  cannot  be  separated  from  the  clutter  unless 
either  its  cross-section  is  much  greater  than  20  dBsm  or  its 
Doppler  is  significantly  different  from  that  of  the  clutter.  The 
limitations  of  Doppler  filtering  to  improve  S/C  were  indicated  in 
Section  V-A. 

An  aircraft  cross-section  is  on  the  order  of  15  dBsm  so 
that  the  maximum  S/C  for  a stationary  aircraft  is  on  the  order  of 
only  -5dB.  The  S/C  will  be  much  less  with  practical  beamwidths 
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which  contain  clutter  areas  with  dimensions  large  compared  to  those 
of  the  aircraft.  Consequently,  detecting  moving  aircraft  will  be 
possible  only  if  their  Doppler  is  separated  sufficiently  from  that 
of  the  clutter. 

5.4  Sample  Waveforms  from  Range-Doppler  Maps 

Sample  waveform  designs  based  on  the  R-R  maps  of  the 
clutter  and  aircraft  for  a few  bistatic  configurations  are  now 
considered.  It  is  noted  that  resolution  and  ambiguity  for  a 
bistatic  and  monostatic  radar  are  defined  in  terms  of  PRF,  band- 
width and  integration  time  in  the  same  way. 

5.4.1  Low  Altitude  Satellites 

The  first  example  is  for  a bistatic  radar  on  a pair 
of  co-orbital  low  altitude  satellites  whose  orbit  plane  is  assumed 
adjusted  to  permit  surveillance  of  suitable  portions  of  the  fence. 
(Several  pairs  of  satellites  would  be  needed  for  continuous  cover- 
age because  of  the  limited  horizon  of  low  altitude  satellites.) 

Figure  V-2  illustrates  the  low  altitude  bistatic 
radar  for  this  example.  The  antenna  on  each  satellite  is  assumed 
to  have  a pencil  beam  and  the  two  beams  are  equal.  Since  satellite 
S,  is  closest  to  the  view  point  x y , its  beam  intersection  with 
the  ground  will  be  entirely  within  that  of  satellite  S2  and  hence 
will  define  the  illuminated  clutter  region  as  indicated  in  Fig.  V-3 
The  range  and  Doppler  distribution  of  the  clutter  within  this  regirn 
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Fig.  V-2.  Low  altitude  bistatic  radar. 


nical 
I exagge 


can  be  determined  from  the  mapping  of  the  clutter  boundary  onto 
the  R-R  plane  as  shown  in  Fig.  V-4  for  a beainwidth  of  0.5°.  The 
contour  shown  is  assumed  to  be  along  the  first  null  of  the  beam. 
It  is  also  assumed  that  the  beam  sidelobes  are  sufficiently  deep 
that  any  clutter  contribution  through  them  can  be  neglected. 


Also  shown  in  Fig.  V-4  is  the  R-R  mapping  of 
0.3  km/s  (~  Mach  .9)  aircraft  on  the  beam  boresight  with  headings 
relative  to  the  satellites'  velocity  (see  Fig.  V-2)  of  between  0 
and  360  degrees.  Each  relative  heading  corresponds  to  a particular 
Doppler  along  the  line  indicated  in  Fig.  V-4. 


The  relative  clutter  density  in  the  R-R  plane  is 
given  in  Fig.  V-5  for  arbitrary  points  x,  y on  the  earth.  For  the 

given  viewing  geometry  with  x - y = 500  km,  J = 160  seconds.  J, 

(j)  (p 

is  approximately  constant  within  illuminated  regions  a few  tens  of 
kilometers  in  extent  as  indicated  by  the  5°  beam  interse'*  :>n  in 
Fig.  V-5.  For  large  illuminated  regions  will  vary  over  the 
illuminated  region,  so  it  would  be  appropriate  to  determine 

separately  for  each  range  slice  or  group  of  range  slices.  In  this 


case,  with  a 0.5  beam,  it  is  not  necessary  to  do  this  and  J can 


be  assumed  constant  within  the  beam. 


If  a resolution  cell  is  fully  occupied  by  clutter 
then  selecting  the  viewing  geometry  to  minimize  J,  minimizes  the 
clutter  competing  with  an  aircraft  in  that  cell,  but  may  be  futile 
for  practical  S/C  ratios  as  indicated  in  Section  V-1.  On  the  other 
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Fig.  V-4.  Clutter  and  aircraft  mapped  onto  the  range-doppler 
plane  for  low  altitude  satellites. 


I 


I 


RANGE  (km) 


hand,  recalling  that  0^  = J^dRdR,  we  note  that  with  a fixed  range 
resolution  dR  and  a fixed  beam  size  with  its  corresponding  fixed 
total  clutter  within  the  beam  (except  for  geometry  dependence  of  y) , 
selecting  large  results  in  the  clutter  being  concentrated  in  a 
narrow  Doppler  band.  A limiting  case  is  when  becomes  infinite 
corresponding  to  a zero  clutter  Doppler  and  clutter  Doppler  spread. 
When  the  selected  geometry  also  results  in  the  aircraft  and  clutter 
being  separated  in  Doppler,  then  the  narrow  clutter  Doppler  band 
can  facilitate  clutter  filtering  and  target  detection. 

For  determination  of  waveform  parameters  a 0.5° 
beam  is  assumed.  The  characteristics  of  the  clutter  in  the  beam 
are  summarized  in  Table  V-1. 


The  range  and  Doppler  ambiguities  are  related  by 


= 75 


= § T 

2 p 


(V-1) 


where  R^  = range  ambiguity,  KM 

R = Doppler  ambiguity,  KM/s 

a 

A = wavelength,  M 

T = interpulse  period,  s 

c ,=  velocity  of  light,  3 x 10^  KM/s 
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TABLE  V-1 


CLUTTER  CHARACTERISTICS  FOR 
LOW  ALTITUDE  BISTATIC  RADAR  EXAMPLE 


Doppler  spread 
Range  spread 
Illuminated  area 
Bistatic  incidence  angle 
Jacobian 
J cos  0 

Total  clutter  cross-section 

Clutter  cross-section  in 
clutter  filled  R-R  cell 


AR  = 81  m/s 

s 

AR  = 11  km 

s 

A = 170  km^ 

())  = 51° 

J = 250  sec 

Jg  = 160  sec 

2 

= Y 10  m 

• • 

a = YJodRdR  = 160  y dRdR  m 

c b 


since  the  aircraft  must  be  in  the  beam  for  detection, 
it  will  be  within  the  range  interval  of  the  clutter  as  defined  by 
the  beam  (for  most  practical  geometries) . In  order  to  prevent 
clutter  foldover  or  ambiguous  range  measurements  and  to  provide  as 
large  an  unambiguous  Doppler  interval  as  possible,  the  range  ambiguity 
is  set  equal  to  the  range  extent  of  the  clutter  in  the  beam.  The 
inter-pulse  period  and  the  Doppler  ambiguity  are  found  from  Eq.  V-1. 
The  basic  waveform  parameters  are  summarized  in  Table  V-2  for 
Y = 0.03,  0.1  and  0.3  ra . 

Since  the  clutter  density  is  approximately  constant 
within  the  beam,  that  range  cell  which  has  the  greatest  spread  in 
clutter  Doppler  will  have  the  greatest  clutter  cross-section.  For 
the  example  under  consideration  this  is  approximately  the  range  cell 
though  the  beam  boresight  (see  Fig.  V-4)  which  has  an  R spread  of 
approximately  78  m/s.  The  clutter  within  that  range  cell  is,  for  a 
fairly  conservative  value  of  y = 0.1, 

Lo^  = 3.7  10^  m^  = 45  dBsm. 

These  clutter  characteristics  are  inclided  in  Table  V-2.  Other 
range  cells  have  smaller  Doppler  spreads  and  hence  have  snialler 
clutter  cross-sections  and  greater  clear  areas. 

When  the  clutter  Doppler  spread  is  a small  fraction 
of  the  Doppler  ambiguity  spacing  there  is  a substantial  clear  aroa 
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TABLE  V-2 


WAVEFORM  AND  CLUTTER  PARAMETERS 
FOR  LOW  ALTITUDE  SATELLITE  EXAMPLE 


Waveform 


Wavelength,  X , m 

0.3 

0.1 

0.03 

Range  ambiguity,  R^  , km 

11.0 

11. n 

11.0 

Pulse  period,  Tp  , ys 

74 

74 

74 

• 

Doppler  ambiguity,  R^  , m/s 

2000 

680 

200 

Pulse  repetition  rate,  PRF  , 

KHz 

13.5 

13.5 

13.5 

Range  resolution,  AR,  m 

30 

30 

30 

Pulse  length,  x , ys 

0.2 

0.2 

0.2 

Doppler  spread/Ambiguity  spacing. 

% 

4 

12 

39 

Clear  area  % 

96 

88 

61 

Clutter 

Y ^ 0.1 

Total  Clutter  ~ 

10^ 

2 

m 

= 70. 

4 dBsm 

Max.  Doppler  velocity 
spread  in  range  cell  ^^SC 

= 

78 

m/s 

Max.  clutter  in 

range  cell  Aa 

-- 

3. 

7 X 10^ 

m =46 

dBsm 
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RELATIVE  DOPPLER  R,^^^  I km/sec) 


18-3-19111 


Fig.  V-6.  Target-clutter  relative  doppler  for  low  altitude 
satellite  example. 
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in  the  R-R  plane  free  of  clutter.  This  fraction  and  the  correspond- 
ing fraction  of  clear  area  are  also  included  in  Table  V-2  as  percent- 
ages. A longer  wavelength  results  in  a greater  percentage  of  clear 
area,  however,  a correspondingly  larger  aperture  is  then  required 
to  keep  the  beamwidth  constant. 

Only  the  example  with  y = 0.1  is  subsequently 
considered.  For  the  aircraft  to  be  separate  from  the  clutter  in 
the  central  range  cell,  the  target  R must  differ  from  that  of  the 
clutter  by  at  least  half  of  the  clutter  R spread,  or  39  m/s.  Figure 
V-6  shows  the  difference  in  R between  an  aircraft  on  the  beam 
boresight  and  the  co-located  clutter  as  a function  of  the  relative 
aircraft  heading  (See  Fig.  V-2).  The  maximum  relative  R is 
210  m/s,  well  within  the  Dopp'er  ambiguities.  It  is  evident  that 
the  relative  R will  exceed  39  m/s  for  relative  headings  from  -30® 

to  135®  or  for  about  90%  of  all  relative  headings  of  interest  (i.e., 

• • 

aircraft  with  only  positive  R or  by  symmetry  with  only  negative  R) . 

(Relative  headings  180°  from  these  also  satisfy  ^pler  separa- 

tion constraint  but  in  this  example  they  would  represent  aircraft 
traveling  away  from  the  U.  S.  and  henc«  pose  no  threat.)  In  either 
case  the  magnitudes  of  the  relativ  . R s are  the  same.  For  aircraft 
at  other  locations  in  this  range  cell  the  relative  Doppler  would  be 
little  different  because  the  dimensions  of  the  beam  intersection 
with  the  earth  are  small  compared  to  the  ranges  to  the  satellites. 

A reduction  in  the  allowable  range  of  relative  headings  due  to 
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waveform  weighting  will  be  discussed  below.  The  orientation  of  the 
satellite  orbit  relative  to  potential  aircraft  compass  headings 
should  be  selected,  if  possible,  so  that  the  aircraft  will  always 
be  separate  from  the  clutter  in  Doppler. 

The  relation  of  the  clutter,  aircraft  and  waveform 
ambiguities  for  an  aircraft  on  the  beam  boresight  with  maximum 
relative  Doppler  is  summarized  in  Fig.  V-7. 

To  obtain  a 15  dB  signal/clutter  ratio  in  this 
example,  45  dB  of  clutter  rejection  is  required.  In  order  to  obtain 

the  necessary  Doppler  sidelobes  some  form  of  waveform  weighting  is 
required.  However.,  since  the  beam  limits  the  range  extent  of  the 

clutter  and  the  PRF  is  selected  to  prevent  clutter  foldover  in 
range,  only  the  waveform  response  at  zero  relative  range  need  be 
considered.  (See  Fig.  V-7).  One-way  Hamming  weighting  can  provide 
-42  dB  sidelobes  with  only  1.34  dB  of  mismatch  loss.  Figure  V-8 
shows  a cut  along  the  Doppler  axis  of  the  ambiguity  function  of  a 
32  pulse  one  way  Hamming  weighted  burst.  Except  for  the  first  four 
sidelobes  they  art;  all  below  the  desired  45  dB.  The  clutter  is 
spread  over  1560  Hz  which  corresponds  approximately  to  the  width 
of  the  region  containing  the  four  sidelobes  exceeding  -45  dB.  The 
worst  effect  resulting  from  clutter  located  near  the  target  or 
ambiguity  as  indicated  in  Fig.  v-7  would  be  to  reduce  the  S/C  by 
less  than  about  3 B.  On  this  basis  the  sidelobes  can  reject  the 
clutter  when  it  is  between  1 and  12.5  KHz  (52  and  630  m/s).  Since 
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Fig.  V-7.  Ambiguity  locations  for  low  altitude  satellite 
example. 
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the  clutter  Doppler  spread  is  78  m/s  an  aircraft  clutter  separation 
greater  than  52  + 78/2  = 91  m/s  will  assure  adequate  clutter 
suppression.  From  Fig.  V-6  it  is  seen  that  a relative  Doppler 
separation  of  91  m/s  corresponds  to  target  headings  between  -11° 
and  117°,  or  about  70%  of  the  possible  headings  of  interest. 

The  requirement  that  the  peak  sidelobes  of  the 
ambiguity  function  be  below  -45  dB  is  conservative  since  the  clutter 
Doppler  spread  also  includes  sidelobe  nulJs.  The  integrated  response 
more  accurately  characterizes  the  response  of  the  waveform  to  clutter. 
Figure  V-9  gives  the  integrated  response  of  the  above  waveform  for 
various  locations  of  the  center  of  the  clutter  Doppler  spread. 

When  the  clutter  is  at  zero  frequency  the  integrated  clutter  is 
nearly  10  dB  below  that  if  the  ambiguity  response  were  unity 
throughout  the  region  occupied  by  the  clutter.  The  net  effect  is 
that  the  integrated  response  is  significantly  less  than  that 
indicated  by  the  sidelobes  on  the  Doppler  axis.  Since  the 
integrated  response  is  less  than  -50  dB  between  ambiguities,  there 
will  be  a 5 dB  margin  of  clutter  rejection  with  this  waveform. 

The  duty  cycle  of  the  waveform  can  be  increased 
(for  higher  average  power)  without  altering  the  ambiguity  response 
on  the  Doppler  axis  by  using  an  FM  pulse  with  a 5 MHz  bandwidth, 
the  same  bandwidth  as  the  original  unmodulated  0.2  ys  pulse. 

The  choice  of  32  pulse  Hamming  weighting  requires  an 
integration  time  of  32  j PRF  - 2.4  ms.  The  Doppler  resolution  is 
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20  m/s  and  the  corresponding  coverage  rate  is 


= 170  km^ 

dt'  2.4  ms 


7 X 10^  km^/s 


3 2 

which  exceeds  the  2.9  x 10  km  /s  rate  required  to  search  the  entire 
fence  in  one  hour  (Section  III) . Several  tradeoffs  are  possible 
to  reduce  the  coverage  rate  to  the  minimum  that  is  necessary  and 
thereby  use  the  minimum  power-aperture  as  indicated  in  Fig.  III-l. 

5.4.2  High  Altitude  Satellites 

A second  example  of  the  application  of  R-R  maps  to 
waveform  design  is  for  the  bistatic  radar  with  the  satellites' 
angular  rates  at  apogee  equal  to  that  of  the  earth.  The  coverage 
by  this  radar  was  given  in  Fig.  IV- 3.  Since  the  satellites'  rates 
at  apogee  match  that  of  the  earth,  the  clutter  Doppler  is  zero  as 
with  a geostationary  satellite.  The  principal  differences  between 
this  satellite  configuration  and  that  of  a geostationary  one  are 
that  the  bistatic  geometry  is  significantly  different  providing 
different  bistatic  aspect  angles  and  that  the  12  hour  period 
satellites  are  not  on  station  continuously.  Also  as  a result  of 
the  satellites'  motions  the  clutter  is  not  zero  at  all  times. 

R-R  clutter  maps  of  the  surveillance  fence  shown  in 
Fig.  V-10  with  corresponding  R-R  maps  of  aircraft  at  various  compass 
headings  within  the  fence  are  presented  in  Figs.  A-4  to  A-11  of 
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Appendix  A.  The  surveillance  fence  is  divided  into  a grid  5 
longitude  by  1°  latitude  and  each  grid  point  mapped  into  the  R-R 
plane  to  characterize  the  clutter.  Aircraft  at  each  grid  point 
with  an  altitude  of  10  km  and  moving  horizontally  at  0.3  km/s  on 
compass  headings  of  135,  180  and  225  degrees  are  also  mapped  into 
the  R-R  plane.  These  headings  bound  those  associated  with  the 
representative  routes  described  in  Section  I.  The  Doppler  spread 
of  the  clutter  or  the  aircraft,  over  the  extent  of  the  fence,  is 
small  (<15  m/s)  so  that  their  relation  in  the  R-R  plane  is 
approximately  that  of  the  line  segments  shown  in  Fig.  V-11.  For 
Fig.  V-11  the  satellites  are  at  apogee  and  1 and  3 hours  past  apogee. 
The  clutter  and  aircraft  are  seen  to  be  widely  separated  in  Doppler 
(>70  m/s)  for  periods  of  at  least  6 hours  (3  hours  before  and  3 
hours  after  apogee) . 

The  intersection  of  a 0.1®  beam  and  the  earth's 
surface  is  mapped  in  Fig.  V-12  for  point  I of  the  fence  for  1 hour 
and  3 hours  after  apogee.  (See  Fig.  V-10) . At  apogee  the  satellite 
rates  match  the  earth  rate  so  that  the  clutter  Doppler  is  zero  and 
its  spread  is  zero.  The  Doppler  spread  across  the  beam  and  hence 
the  clutter  Doppler  spread  in  a range  cell  is  seen  to  be  very  small 
(<0.6  m/s)  for  this  example.  The  total  clutter  cross-section  and 
cross-section  in  a bistatic  range  cell  can  be  found  from  maps  like 
this  and  the  clutter  density.  (However,  the  actual  clutter  Doppler 
spread  may  be  on  the  order  of  a few  meters  per  second  due  to  motion 
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of  individual  scatterers  comprising  the  clutter.  The  rms  clutter 
velocity  for  a monostatic  radar  is  typically  < 4 m/s^^^  so  that 
for  a bistatic  system  it  can  be  expected  to  be  in  the  order  of 
Vcb  1 4 cos  D/2  m/s  where  B is  the  bistatic  angle.)  The  bistatic 
range  extent  of  the  clutter  is  used  to  determine  the  waveform  period 
as  was  done  with  the  lew-altitude  satellite  example. 

The  variation  of  the  clutter  density  over  the  fence 
is  illustrated  in  Figs.  A-12  and  A-13  of  Appendix  A for  this  example. 
The  clutter  density  and  clutter  cross-section  in  a range  slice  for 
tne  beam  boresighted  on  point  I are  tabulated  in  Table  V-3  which 
summarizes  the  clutter  characteristics  and  possible  waveform 
parameters.  Possible  techniques  for  obtaining  large  clutter  re- 
jection are  discussed  in  Section  Vi  .M+-hough  the  clutter  distri- 
bution (Fig.  V-11)  varies  with  time,  the  clutter  range  extent, 
which  determines  the  waveform  characteristics,  does  not  vary  greatly 
throughout  the  fence  nor  with  time  so  that  a single  waveform  could 
be  applicable  independent  of  time.  However,  both  the  range  and 
Doppler  location  of  the  filters  would  need  to  be  shifted  as  a 
function  of  time. 

The  maximum  dwell  time  in  each  beam  position,  and 
hence  the  maximum  number  of  pulses  N that  can  be  integrated 
consistent  witn  the  reauired  coverage  rate  is 

„ ^ illumination  area 

N < P RF 

coverage  rate 
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and  is  included  in  Table  V-3.  N is  very  large  and  will  not  limit 
attainment  of  the  required  clutter  filter  characteristics. 

The  aircraft  cross-section  is  aependent  on  the 

bistatic  viewing  geometry  and  requires  four  angles  for  its 

description.  These  four  angles  are  the  bistatic  angle,  the  bistatic 

pitch,  the  bistatic  aspect  and  the  bistatic  roll  which  are  illustrated 

in  Appendix  A,  Figs.  A-2  and  A-3.  The  variations  of  these  angles 
over  the  various  aircraft  locations  at  the  5°  x 1°  grid  points 

within  the  surveillance  fence  are  tabulated  in  Table  V-4.  Three 
aircraft  headings  at  each  grid  point  are  considered;  135° (south- 
east) , 180°  (due  south)  and  225°  (south-west)  . The  avera-jS, 
maximum,  minimum  and  spread  (difference  between  the  iraximum  and 
minimum)  of  the  four  angles  is  given  for  the  jatelliteK-.  at  apogee 
(T  = 0 HR)  and  three  hours  later  (T  = 3 HR) . 

Table  V-4  indicates,  for  this  example,  that  although 
the  bistatic  angle  for  the  viewing  of  aircrait  anywhere  within  the 
surveillance  region  is  nearly  constant  (and  not  very  large) , the 
bistatic  pitch,  aspect  and  roll  have  significant  variations.  Thus, 
it  is  expected  that  the  aircraft  cross-section  would  vary  signifi- 
cantly depending  on  its  location  and  heading  and  the  time  of 
observation  or  locations  of  tne  satellites.  Aircraft  maneuvers 
away  from  the  straight  and  level  flight  would  result  in  further 
variations  in  the  bistatic  viewing  angles.  The  effect  on  aircraft 
bistatic  viewing  angles  of  the  bistatic  radar  on  satellites  on  other 
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TABLE  V-3 


WAVEFORM  AND  CLUTTER  PAIAMETERS  FOR 
HIGH  ALTITUDE  SATELLITES 


Time  from  apogee  - hr 

0 

1 

3 

CLUTTER 

R extent,  m/s 

0.14 

0.6 

Range  extent,  km 

35 

34 

25 

• 

F-R  Clutter  density,  dbs 

00 

66 

58 

* ( cl ) 

R spread,  m/s 

~0 

0.03 

0.2 

Clutter  0,  dbsm^^^ 

56 

56 

56 

Clutter  rejection 

56 

56 

56 

Bistatic  angle,  deg 

60 

60 

80 

WAVEFORM 

Period,  ys 

230 

prf,  Hz 

4400 

• 

R ambiguity,  m/s 

215 

Max.  no.  integrated  pulses 

(d) 

5700 

(a)  In  30  m range  slice 

(b)  Required  for  15  db  S/C  on  15  dbsm  target 

(c)  Range  ambiguity  equal  to  clutter  extent;  X = 0.1  m 

(d)  From  maximum  dwell  time  consistent  with  ground  coverage  rate 
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TABLE  V-4 


SUMMARY  OF  BISTATIC  AIRCRAFT  VIEWING  ANGLES* 


Time  After  Apogee 

T = 

0 HR 

T = 

3 HR 

Heading 

= 135 

o 

Angles,  Deg. 

Average 

Spread 

Max 

Min 

Average 

Spread 

Max 

Min 

Bistatic  Angle 

58 

1 

59 

58 

78 

2 

79 

77 

" Pitch 

39 

80 

70 

-10 

39 

77 

69 

-9 

" Aspect 

195 

111 

250 

139 

205 

112 

261 

149 

" Roll 

74 

148 

125 

-23 

73 

149 

126 

-24 

Heading 

= 180 

o 

Bistatic  Angle 

58 

1 

59 

58 

78 

2 

78 

77 

" Pitch 

0 

98 

-49 

49 

1 

96 

49 

-47 

" Aspect 

160 

31 

178 

148 

174 

30 

189 

159 

" Roll 

90 

97 

139 

41 

89 

100 

138 

38 

Heading 

= 225' 

0 

Bistatic  Angle 

58 

1 

59 

58 

78 

2 

79 

77 

" Pitch 

-39 

80 

10 

-70 

-38 

80 

11 

-69 

" Aspect 

195 

111 

250 

139 

203 

107 

256 

149 

" Roll 

129 

148 

203 

55 

126 

146 

199 

53 

*For  the  satellite  orbits  of  the  high  altitude  example  of 
Section  V-5.3  and  Figure  IV-3.  See  also  Appendix  A, 
Figures  A-2  and  A-3. 
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orbits  needs  investigation.  However,  it  is  expected  that  if  a large 
fence  is  to  be  surveyed,  then  the  changes  in  the  aircraft  bistatic 
viewing  angles  will  also  be  large.  It  is  also  expected  that  the 
variation  of  bistatic  angle  for  lower  altitude  satellites  would  be 
greater  than  for  the  example  presented.  It  is  also  noted  that 
bistatic  measurements  of  aircraft  cross~section  require  combinations 

of  four  different  angles  which  involves  an  extensive  number  of 
observations.  For  example,  if  only  five  different  values  of  the 

four  bistatic  angles  were  selected  (which  is  very  unlikely  to  be 

4 

adequate)  then  5 or  625  separate  measurements  are  required. 

The  sample  waveform  designs  of  this  section  have 
served  to  illustrate  some  of  the  characteristics  and  applications 

of  bistatic  R-R  maps  of  clutter  and  aircraft  and  of  clutter  density 

* 

in  the  R-R  plane.  They  demonstrate  the  possibility  of  selecting 
relatively  simple  waveforms  that  could  provide  aircraft  detection 
consistent  with  the  clutter  filtering.  The  possible  need  for 
tracking  the  center  of  the  filters  in  both  range  and  Doppler  v^as 
illustrated.  For  both  the  low  and  high-altitude  satellite  examples 
uhe  target  and  clutter  were  separated  so  that  detection  is  possible 
provided  Doppler  filters  have  sufficiently  low  sidelobes  for  the 
assumed  bistatic  cross-sections.  There  is  a need  for  wideband  wave- 
forms at  large  bistatic  angles  if  there  is  to  be  maximum  clutter 
rejection.  Based  on  these  two  examples,  it  is  not  expected  that 
waveforms  will  be  a limiting  factor  in  the  design  of  satellite-borne 
bistatic  radar  for  detection  of  aircraft.  Clutter  filters  are 
considered  in  the  next  section. 
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FILTER  DESIGN 


VI . 

6.1  Waveform  Weighting 

Waveform  weighting  is  employed  to  redu::e  range  and  Doppler 
sidelobes  to  levels  sufficiently  low  that  clutter  can  be  rejected. 

An  application  of  Hamming  weighting  is  given  in  Section  V.  An 
example  of  viewing  the  earth  from  40,000  km  that  resulted  in  a 
clutter  cross-section  of  ~60  dBsm  and  -60  dB  sidelobe  level  is 
given  in  Section  III.  Since  the  need  of  sidelobe  levels  below 
that  attainable  with  Hamming  weighting  may  arise,  other  weightings 
will  be  briefly  considered  here. 

Weighting  results  in  a signal/noise  ratio  loss  when  compared 
to  an  unweighted  waveform.  Since  the  signal  is  coherently  integrated, 
the  signal  energy  is  proportional  to  the  square  of  the  sum  of  weights 
and  the  noncoherent  noise  energy  is  proportional  to  the  sum  of  the 
squares  of  the  weights.  Therefore,  the  S/N  for  a weighted  waveform 
(on  receive  only)  is  proportional  to 


n=l 
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where  W(n)  are  the  weights  and  N the  number  of  pulses.  The  S/N  for 
the  unweighted  waveform,  W(n)  =1,  is  proportional  to  the  number  of 
pulses  N,  so  that  the  S/N  losses  due  to  one-way  weighting  are 


For  unequal  weighting  on  both  transmit  and  receive  the  loss  is 


This  loss  will  be  determined  when  considering  various  weightings. 
6.1.1  Hamming  Weighting 

Hamming  weighting  on  receive  only  results  in  -42  dB 
sidelobes.  The  sidelobe  level  is  insensitive  to  the  number  of 
pulses  over  which  the  weighting  is  made  as  long  as  the  number  of 
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pulses  provides  an  adequate  approximation  to  the  integral  of  the 
Hamming  weighting  function.  Figure  V-8  shows  the  response  for 
weighting  over  32  pulses.  The  S/N  loss  is  only  1.34  dB. 

Lower  sidelobes  result  with  weighting  on  both 
transmit  and  receive.  Figure  VI-1  is  the  Doppler  axis  cut  of  the 
ambiguity  function  for  two-way  Hamming  weighting  over  32  pulses. 

The  sidelobe  levels  are  about  20  dB  lower  than  for  one-way  weighting 
and  the  loss  is  only  2.6  dB. 

6.1.2  Weights  from  Digital  Filter  Design 

It  is  possible  to  design  digital  finite  impulse 

(T) 

response  filters  with  very  low  sidelobes.  Using  such  a filter 
to  process  the  received  waveform  can  be  expected  to  result  in 

/ 0 \ 

substantial  clutter  reduction  capability.  Using  a computer  program' 
for  the  design  of  optimal  finite  impulse  response  (FIR)  linear  phase 
filters,  some  filter  responses  have  been  calculated.  One  example 
is  given  in  Fig.  vi-2.  The  filter  specifications  and  resulting 
weights  are  given  in  Table  VI-1.  The  depth  of  the  sidelobes  is  a 
functiwn  of  filter  length  (the  number  of  pulse  weighted) , and  the 
widths  of  the  pass,  stop  and  transition  bands.  A wide  transition 
band  tends  to  result  in  smaller  sidelobes  but  will  also  tend  to 
reduce  the  clear  area  between  ambiguities.  As  indicated  in  Section 
V a small  clear  region  restricts  the  aircraft  headings  for  which  the 
clutter  can  be  rejectee. 
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TABLE  VI-1 


PULSE  WEIGHTING  DERIVED  FROM  DIGITAL  FILTER  DESIGN: 

-55db  Sidelobes 


FINITE  IMPULSE  RESPONSE  (FIR) 
LINEAR  PHASE  DIGITAL  FILTER  DESIGN 
REMEZ  EXCHANGE  ALGORITHM 


BANDPASS  FILTER 


FILTER  LENGTH  = 32 


*****  IMPULSE  RESPONSE  ***** 


H( 

1) 

= 

-0.21792529E-02 

H( 

32) 

H( 

2) 

-0.33740327E-02 

=: 

H( 

31) 

H ( 

3) 

= 

-0.51930882E-02 

z= 

H( 

30) 

H( 

4) 

= 

-0.67545921E-02 

= 

H( 

29) 

H{ 

5) 

= 

-0.73748156E-02 

= 

H( 

28) 

H( 

6) 

zs 

-0.62420666E-02 

=: 

H( 

27) 

H( 

7) 

5= 

-0.25451574E-02 

=5 

H( 

26) 

H( 

8) 

= 

0.43474101E-02 

=s 

H( 

25) 

H( 

9) 

zz 

0.14719278E-01 

= 

K( 

24) 

H(10) 

= 

0.28370515E-01 

H( 

23) 

H(ll) 

= 

0.44553712E-01 

= 

H( 

22) 

H(12) 

= 

0.62006533E-01 

= 

H( 

21) 

H(13) 

=5 

0.79088032E-01 

H( 

20) 

H(14) 

= 

0.93998909E-01 

= 

H( 

19) 

H(15) 

0.10505038E+00 

=r 

H( 

18) 

H(16) 

= 

0.11093324E+00 

H( 

17) 

LOWER  BAND  EDGE 
UPPER  BAND  EDGE 
DESIRED  VALUE 
WEIGHTING 
DEVIATION 
DEVIATION  IN  DE 


BAND  1 
0,0 

0.020000000 

1.00000000 

1.00000000 

0.018813349 

0.161890686 


BAND  2 
0.100000024 
0.500000000 
0.0 

10.0000000 

0.001881335 

-54.5106659 
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The  ambiguity  response  along  the  Doppler  axis  for  a 
receiver  using  this  filter  is  shown  in  Fig.  Vl-i.  Its  sidelobe 
response  (-55  dB)  is  signif icantly  lower  than  that  for  Hamming 
weighting  but  at  the  expense  of  4.5  dB  S/N  loss.  This  loss  is  due 
in  part  to  the  fact  that  negative  weights  are  used  to  obtain  the 
low  sidelobes. 

The  integrated  response  of  this  ambiguity  function 
to  the  clutter  distribution  associated  with  the  low  altitude 
satellite  example  of  Section  V is  shown  in  Fig.  VI-4.  The  clutter 
rejection  for  this  distribution  of  clutter  is  ~ 60  dB. 

\nother  example  of  an  FIR  filter  design  and  associated 
ambiguity  function  is  given  in  Table  VI-2  and  Fig.  VI-5,  The  side- 
lobes  are  less  than  -70  dB  but  the  S/N  loss  is  7 dB. 

6.2  Implementation  Considerations 

The  sample  filter  designs  discussed  above  indicate  that 
substantial  sidelobe  reduction  and  hence  clutter  rejection  is 
possible  but  only  with  significant  S/N  loss.  However,  there  are 
several  important  considerations  with  regard  to  the  implementation 
of  these  designs.  The  filter  design  should  be  optimized  for  the 
specific  application  both  with  regard  to  integrated  clutter  response 
and  S/N  loss.  Digital  signal  processing  logic  needs  to  be  formulated 
to  efficiently  provide  the  desired  filtering  and  target  range  and 
Doppler.  The  stability  and  accuracy  requirements  of  receiver 
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RESPONSE  (dB) 


FREQUENCY  (Hz) 


Ficf.  VI-3.  Ambiguity  function  on  doppler  axis  for  weighting 
of' Table  VI-1. 
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TABLE  VI- 2 


PULSE  WEIGHTING  DERIVED  FROM  DIGITAL  FILTER  DESIGN: 

-70db  Sidelobes 


FINITE  IMPULSE  RESPONSE  (FIR) 
LINEAR  PHASE  DIGITAL  FILTER  DESIGN 
REMEZ  EXCHANGE  ALGORITHM 


BANDPASS  FILTER 

FILTER  LENGTH  = 32 

*****  impulse  RESPONSE  ***** 


H( 

1) 

0.60465327E-03 

= 

H( 

32) 

H( 

2) 

= 

0.14172327E-02 

= 

H( 

31) 

H( 

3) 

= 

0. 22057937E-02 

= 

H( 

30) 

H( 

4) 

= 

0.20479802E-02 

= 

H( 

29) 

H( 

5) 

= 

-0.16512611E-03 

= 

H( 

28) 

H( 

6) 

-0.51321425E-02 

H( 

27) 

H( 

7) 

= 

-0.12332562E-01 

= 

H( 

26) 

H( 

8) 

s: 

-0.19445047E-01 

=s 

H( 

25) 

H( 

9) 

s 

-0.22437945E-01 

= 

H( 

24) 

H(10) 

= 

-0.16578048E-01 

= 

H( 

23) 

H(ll) 

s 

0.18148581E-02 

St 

H( 

22) 

H(12) 

a 

0.33475608E-01 

s: 

H( 

21) 

H(13) 

= 

0.75001478E-01 

H( 

20) 

H(14) 

= 

0.11904991E+00 

H( 

19) 

H(15) 

= 

0.15610170E+00 

= 

H( 

18) 

H(16) 

= 

0.17729408E+00 

H( 

17) 

BAND  1 

LOWER  BAND  EDGE  0.0 

UPPER  BAND  EDGE  0.050000001 

DESIRED  VALUE  1.00000000 

WEIGHTING  1.00000000 

DEVIATION  0.014156718 

DEVIATION  IN  DB  0.122098148 


BAND  2 
0.149999976 
0.500000000 
0.0 

50.0000000 

0.000283134 

-70.9601593 
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-12,000  -8000  -4000  0 4000  8000  12,000 


FREQUENCY  (Hz) 

Pig,  VI-5.  Ambiguity  function  on  doppler  axis  for  weighting 
of  Table  VI-2. 
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hardware  preceding  the  digital  processing  may  be  high.  These 
requirements  need  to  be  defined  and  needed  hardwire  developments 
identified.  The  possibility  of  accurately  controlled  transmitter 
weighting  should  not  be  overlooked.  The  transmit  modules  proposed 
for  constructing  large  planar  arrays  may  accommodate  accurate  level 
setting.  The  high  output  solid  state  devices  used  in  these  modules 
contain  many  parallel  sections  and  the  switching  of  these  may  be 
one  way  to  permit  accurate,  level  setting. 
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VTI. 


SUMMARY  AND  RECOMMENDATIONS 


Satellite-borne  bistatic  radars  can  be  placed  in  orbits 
so  that  flying  aircraft  and  ground  clutter  do  not  overlap  in 
Doppler.  Simple  waveforms  can  then  be  designed  which  permit 
detection  of  the  aircraft  using  Doppler  filtering  techniqutj. 
Mappings  of  aircraft  targets  and  ground  clutter  onto  the  .gi  static 
range-Doppler  plane  form  the  basis  of  waveform  design,  Doppler 
filter  capable  of  the  order  of  60  dB  of  clutter  rejectic  •.  is 
required.  Weighting  of  the  waveform  and  filter  to  achir.  a this 
degree  of  clutter  rejection  may  result  in  a loss  of  se'  .-al  dB 
in  signal-to-noise  ratio.  Stringent  stability  and  acc  ;racy 
requirements  are  imposed  on  the  radar  signal  system  nc  <ermit  such 
filters  to  be  used.  Weighting  on  transmit,  using  lev-.'!l  setting 
in  transmit  array  modules,  may  be  useful  as  an  aid  to  obtaining 
required  clutter  rejection. 

Whether  the  ratio  of  aircraft  cross-section  to  clutter 
cross-section  is  larger  for  bistatic  geometry  than  for  monostatic 
would  have  to  be  determined  by  measurements.  For  geometries  with 
large  bistatic  angles,  the  bistatic  radar  must  have  a significantly 
larger  bandwidth  than  the  monostatic  system  for  the  same  ground 
resolution.  For  maximum  clutter  reduction  the  bandwidth  should 
change  with  viewing  geometry  changes. 

The  coverage  characteristics  of  bistatic  radars  on 
satellites  in  representative  orbits  are  described.  For  some  orbits 
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the  bistatic  coverage  is  unique  in  that  it  cannot  be  duplicated 
by  a satellite-borne  monostatic  radar.  For  example,  portions  of 
a coverage  of  a bistatic  radar  can  be  made  equivalent  to  a mono- 
static  radar  fixed  in  inertial  space.  The  possibility  of  bistatic 
range  gating  of  clutter  is  discussed. 

The  variation  of  bistatic  viewing  angles  of  an  aircraft 
crossing  various  portions  of  a fence  may  be  large  so  that  large 
variations  in  cross-section  are  expected.  A complete  description 
of  the  bistatic  cross-section  of  an  aircraft  would  require  extensive 
measurements. 

The  principal  system  requirements  of  a North  American 
surveillance  fence  for  the  detection  of  bombers  are  identified  for 
a space-borne  bistatic  surveillance  and  tracking  rader.  Power- 
aperture  products  consistent  with  a 60  meter  diameter  aperture  and 
1 kw  average  power  are  needed  for  a geostationary  bistatic  system. 
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APPENDIX  A 


Birtcitic  R-R  Maps 

In  this  appendix  the  derivations  of  bistatic  R-R  maps  and 
related  bistatic  parameters  are  outlined.  Both  the  receiver  and 
transmitter  of  the  biscatic  radar  are  assumed  to  be  on  satellites 
which  are  on  elliptic  processing  Kepler ian  orbits.  The  bistatic 
radar  observes  clutter  arising  from  the  surface  of  a rotating 
spherical  earth.  Target  aircraft  moving  above  the  earth's  surface 
may  have  horizontal  and  vertical  velocity  components  relative  to 
the  earth  and  travel  at  arbitrary  headings  or  along  great  circle 
routes.  The  transmitter  and  receiver  beams  may  be  modeled  to  have 
conical  or  elliptical  shapes  and  two  or  four  monopulse  feeds.  A 
highly  flexible  computer  program  has  been  written  incorporating 
these  features  to  obtain  bistatic  R-R  maps  and  other  relevant 
parameters  for  evaluating  bistatic  systems. 

A. 1 Bistatic  Range  and  Range  Gradient 

Figure  A-1  illustrates  a satellite-borne  bistatic  radar 
observing  a point  P and  defines  various  vectors  in  an  inertial 
coordinate  system.  The  vectors  needed  are; 


SI 

satellite 

1 

position 

vector 

S2 

satellite 

2 

position 

vector 

• 

SI 

satellite 

1 

velocity 

vector 

• 

S2 

satellite 

2 

velocity 

vector 
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P - observed  point  position  vector 
P - observed  point  velocity  vector 
R1=P-S1  - range  vector  from  SI  to  P 

R2=P-S2  - range  vector  from  S2  to  P 

Ul=unit  vector  Rl/  |ri| 

U2=unit  vector  R2/  |r2| 

i,  2,  k - tnit  vectors  along  inertial  coordinate 
axis 

The  bistatic  range  is  defined  as 

^ ^ hl|  1^2!  ^ |p-Sl|  + |p-S2| 

B 

2 2 

This  reduces  to  the  monostatic  range  when  the  transmitter  and 
receiver,  i.e.,  the  two  satellites,  are  coincident.  Consequently 
this  derivation  can  also  be  applied  to  monostatic  radars. 

The  bistatic  range  gradient  or  range  direction  is 


For  any  arbitrary  point 

X 

P = y 
z 


= V- 


^(P-Sl)  » (P~S1) 


'(P-S2)  ■ (P~S2) 
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Thus  the  bistatic  range  direction  is  along  the  bisector  of  the 
lines-of-sight  Rl,  R2  to  the  observed  points.  Normalizing  the 
range  gradient  gives  the  unit  vector  UB  along  the  bistatic  range 
direction, 
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UB 


U1  + U2 
|U1  + U2j 


U1  + U2 

[2(1  + U1  • U2)]^/^ 


Defining  the  angle  between  the  line-of-&ight  vectors  as  the 
bistatic  angle  3,  UB  becomes 


UB 


U1  + U2 
2 cos  3/2 


Bistatic  Range  Rate  and  Range  Rate  Gradient 


The  bistatic  range  rate  or  Doppler  is  dR^/dt  or 


• • 


• • 


^B=  t 


[(P  - SI)  • U1  + (P  - S2)  • U2]  . 


The  oistatic  Doppler  gradient  is 


VRb  » V [ 


(P-Sl)  • U1  + (P-S2)  • U2 


3 


i • “ * I!  • >'■2  " ^ + (P-S2)  • Sff  ] 


I [If  -“1  If  • ^ 


I [ 


If.Ul  + II  . U2  + (P-iM  (P-S2). 
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From  the  definitions  of  U1  and  U2,  obtain 


= -L_  fM  - ui  (Ui.  ^)1  . 

|Ri|  [ax  ax] 

Similar  terms  are  obtained  for  the 

partials  of  U2.  For  an  arbitrary  point  P rotating  with  the  earth 


X 

P = y 
z 


and 


where  oi^  is  the  earth  rate.  The  partials  of  P were  given  earlier. 
The  partial  of  P is 


^x 


-0) 


Substitution  into  the  appropriate  terms  above  completes  the  deter- 
mination  of  VRq  at  any  observed  point. 

D 
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A. 2 Bistetic  Angles 

The  cross-section  of  an  object  generally  depends  on  the 
viewing  geometry.  For  a bistatic  radar  there  are  three  parameters 
required  to  describe  the  cross-section  of  an  axi-symmetric  body 
and  four  parameters  for  a body  with  symmetry  in  each  of  two  orthog- 
onal planes  as  is  approximated  by  an  aircraft.  (A  horizontally 
flying  aircraft  has  a horizontal  and  vertical  plane  of  symmetry.) 
The  intersection  of  these  symmetry  planes  defines  the  body  axis. 
Figure  A-2  defines  the  bistatic  angle,  bistatic  aspect  and  bistatic 
pitch  of  an  axi-symmetric  body.  Fox  a body  like  an  aircraft  the 
fourth  parameter  is  the  bistatic  roll.  The  bistatic  roll  is  the 
angle  between  the  bistatic  pitch  plane  and  the  normal  to  the 
"horizontal"  symmetry  plane  of  the  body.  It  is  positive  in  the 
sense  of  a right  hand  screw  along  the  body  axis.  Figure  A- 3 
illustrates  the  bistatic  roll  angle. 

The  bistatic  angles  describing  the  viewing  geometry  of 
areas  (such  as  portions  of  the  earth  surface)  can  be  found  by 
treating  the  surface  as  an  axi-symmetric  body  where  the  surface 
normal  corresponds  to  the  body  axis. 

2^.3  Clutter  Density  in  the  P.-R  Plane 

The  mapping  of  points  on  the  earth's  surface  onto  the 
R-R  plane  has  been  described  above.  However,  clutter  cross-section 
is  area  dependent  so  that  the  clutter  density  in  the  R-R  plane  must 
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Fig.  A-2.  Angles  for  determining  bisbatic  cross  section. 
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be  obtained  by  scaling  elemental  areas.  This  scaling  is  determined 
by  the  Jacobian  of  the  transformation  into  the  R-R  plane. 

The  clutter  cross-section  is  an  elemental  area  given  by 
da  = a*^„dA 

13 

where  a°«  is  the  bistatic  cross-section  density  and  dA  the  elemental 

D 

area  tangent  to  the  earth's  surface.  Expressed  in  terms  of  an 
elemental  area  in  the  R-R  plane  it  is 

^ • 

da  = a _ JdRdR 

o 


where,  for  example,  dR  and  dR  may  be  approximated  by  the  range  and 
Doppler  resolution  of  the  radar.  J is  the  Jacobian  of  the  trans- 

formation  of  dA  to  dR  dR.  Therefore,  the  bistatic  clutter  density 

« 

in  the  R-R  plane  becomes 


da 


BRR 


dR  dR 


= J a 


B 


The  Jacobian  J is  determined  below  from  two  factors;  the  Jacobian 
of  the  transformation  from  x,y  to  latitude,  longitude  coordinates; 
and  the  Jacobian  of  the  transformation  from  latitude,  longitude  to 
R-R  coordinates.  The  first  Jacobian  is  usually  found  by  simpler 
more  direct  means,  but  the  more  general  method  is  used  here  to 
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illustrate  that  method  with  a familiar  case  and  for  consistency 
with  the  determination  of  the  second  factor. 

The  area  dA  can  be  exoressed  as 


dA  = dxdy 


I 


where  x and  y are  related  to  the  latitude  (3 at)  and  longitude  (long) 
of  the  area  dA,  and  the  earth  radius  according  to 

X = lat 

e 


y = Rg  cos  (lat)  long  . 


The  Jacobian  of  this  transformation  is 


*^xyLL 


DET 


dx 

J>JL- 

3lat 


6x 

i^long 

dlong 


9(x,y) 

9(lat,long) 


where  DET  indicates  "determinant  of",  so  that 


■^xyLL  = *^e^ 


^Wl^!| 


and 


dxdy  = ‘^jjyLj  dlat  dlong 


Therefore  in  the  lat,  long  coordinates 


dA  = R cos  (lat)  dlat  dlong  , 
e 


the  familiar  result. 

The  mapping  of  clutter  onto  the  R-R  plane  depends  on  the 
position  and  velocity  of  the  mapped  point,  both  of  which  are 
determined  by  the  latitude  and  longitude  of  the  point.  All  other 
parameters  of  the  transformation  are  independent  of  latitude  and 
longitude.  The  transformation  can  be  represented  by 

R = R (lat,  long) 

• • 

R = R ( ; at , long ) 

where  the  subscript  B,  denoting  a bistatic  system,  is  dropped  for 
compactness  of  notation.  The  Jacobian  of  this  transformation  is 


“^RRLL 


3(R,R) 


3(lat,  long) 
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and  relates  dlat  dlong  and  dR  dR  according  to 


dRdR  = 


dlat  dlong  = J 


RRLL 


LLRR  dlat  dlong  . 


The  Jacobian 


T • 

RRLL 


DET 


dR 

8R 

dlat 


8R 

dlong 

JR 

dlong 


is  evaluated  below. 

Partials  of  R are  found  from 


3H 

, i r?  Ri  + u 

ML] 

3R 

1 

rsjRiL  3|R2  ■\ 

Slat 

2 Ld  lat  3 

lat  J 

' 3 long 

I 

[3 long  3 long J 

8|R1| 
3 lat 

- in 

" • ITit  ' 

dlR2l 

dlat 

= 

n?  3^2 

3|R1| 

dlong 

- Ml  9R1 

■ ' ’dl^ng 

9 

3 1 R^ 
^long 

= 

IT,  3R2 

dlong 

dRl 

01at 


3R2  ^ 3p  ^ „ 
dlat  “ dlat  e 


-sin  lat  cos  (long  + w^t) 
-sin  lat  sin  (long  + uj^t) 
cos  lat 
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9R1 


9 long 


9 R2 
9 long 


9P 


9 long 


= R 


-cos  lat  sin  (long  + ui^t) 
cos  lat  cos  (long  + w^t) 


where  t is  the  time. 


The  partials  of  R are  found  from 


9R 

91at 


J [(P-si) 


9U1 

91at 


9P 

91at 


• U1  + (P-S2)  • 


9U2  . 

9 lat 


9P 

Slat 


9R 

3 long 


3U1  ^ 3P 

3 long  3 long 


U1  + (P-S2) 


3U2 
3 long 


+ 


3P 

3 long 


3U1 

31at 


R1 


3R1 
3 lat 


- R1 


R1 


7 


ilRli 

31at 


3U2 

^lat  ~ 


|R2| 


3R2  , 9R2 

3 lat  “ 9 lat 

R2| 


9U1 
9 long 


|R1| 


3R1 
9 lone 


|R1| 


2 


R1 


ll31i 

3 long 


3U2 
3 long 


|R2| 


3R2 
3 long 


R2 


3R2 
3 long 


R2 


where  the  partials  of  |R1|,  (r2|,  ri  and  R2  are  given  above 
Partials  of  P are  found  from 


P 


^e  ^e 


-sin  (long  + lo^t) 
cos  (long  + w^t) 

0 
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ap 

Slat 


= -0)  R sin  lat 
e e 


-sin  (long  + Wgt) 

cos  (long  + to_t) 

0 

0 


3P 

along 


cos  lat 

e e 


cos  (long  + 

sin  (long  + to^t) 

0 


The  substitution  of  the  partials  determined  above  into 
Jrpll  completes  its  determination.  Combining  *^xyLL 

determines  J and  the  clutter  density  in  the  R-R  plane,  thus 


^ “ '^xyLL^'^RRLL  “ "^xyLL  “^LLRR 


so  that 


BRR  “ ° ^xyLL  ‘^LLRR 


and  the  relative  density  is 

o . 

' _ BRR  _ . 

R ■ gO  " xyLL  LLRR 

An  alternative  relative  R-R  density  results  when  the 
cross-section  density  Yg  relative  to  a projected  area  is  used 
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rather  than  a g.  Thus 


a°g  = Yg  cos  <j) 


where  (j>  is  the  angle  of  incidence  of  the  bistatic  range  vector. 
The  relative  bistatic  density  in  the  R-R  plane  is  then 


BRR 


B 


♦ '^xyLL  ■’lLRR 


and  differs  from  that  initially  derived  by  only  the  cos  (p.  This 
last  form  is  used  in  this  report.  Notice  that  the  dimension  of 
0^  is  seconds. 

A . 4 Example  Mappings  onto  the  R-R  Plane 

To  illustrate  various  R-R  maps  the  bistatic  radar  on  the 
high  altitude  satellites  example  in  Section  V is  used.  These 
satellites  are  on  elliptic,  12  hour  orbits  such  that  their  angular 
rate  at  apogee  matches  that  of  the  earth.  The  clutter  maps  are 
obtained  by  mapping  a 5°  longitude  by  1°  latitude  grid  of  points 

f 

on  the  surveillance  fence  (Fig.  V-10)  into  the  R-R  plane  for  the 
satellites  1 and  3 hours  past  apogee.  Aircraft  at  each  grid  point 
with  an  altitude  of  10  km  and  moving  horizontally  at  0.3  km/s  on 
a heading  of  180  degrees  are  also  mapped  into  the  R-R  plane. 
Because  clutter  or  aircraft  at  different  portions  of  the  fence 
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may  have  the  samr  range  and  Doppler  as  clutter  or  aircraft  in 
other  portions,  there  is  folding  of  the  maps  in  the  R-R  plane. 

To  simplify  the  presentation  of  the  R-R  maps  they  have  been 
divided  into  two  sections  (A  and  B)  in  which  there  is  little  or 
no  folding.  The  clutter  maps  are  given  in  Figs.  A-4,  6,  8 and  10 
and  the  corresponding  aircraft  maps  in  Figs.  A-5,  7,  9 and  11. 

Nine  points  on  the  surveillance  fence  have  been  labeled  as  indi- 
cated  in  Fig.  V-10  to  provide  reference  locations  on  the  R-R  maps. 

Figures  A-12  and  13  show  the  bistatic  relative  cross-section 
density  at  the  grid  points  of  the  fence  when  the  satellites  are  1 
hour  past  apogee.  The  regions  of  large  cross-section  density 
indicate  regions  where  the  Doppler  spread  within  the  beam  will  be 
small.  For  small  beams  where  the  clutter  density  is  constant  within 
the  beam,  the  clutter  cross-section  in  a clutter  filled  R-R  cell  or 
within  the  entire  beam  can  be  directly  found  from  the  density  plots 
and  the  resolution  cell  or  beam  R-R  map. 
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Fig.  A-10.  Clutter  map  of  portion  "B"  of  surveillance  fence: 
t = 3 hr. 
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Fig.  A-12.  Clutter  density  for  portion  "A"  of  surveillance 
fence:  t = 1 hr. 
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Fig.  A-13.  Clutter  density  for  portion  "B"  of  surveillance 
fence;  t = 1 hr. 
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APPENDIX  B 


Range  Gating  Ground  Clutter  in  a Monostatic  Satellite-Borne  Radar 

In  this  appendix  we  investigate  a well  known  technique 
for  eliminating  ground  clutter  in  the  main  lobe  of  a space-borne 
radar  viewing  targets  against  the  background  of  the  earth.  The 
technique  is  to  utilize  range  separation  between  the  target  and 
the  clutter  in  the  beam.  The  range  separation  between  ground 
clutter  and  target  is  investigated  as  a function  of  antenna  beam- 
width,  sensor  scan  angle,  depression  angle,  satellite  altitude, 
target  altitude  and  ground  swath  coverage. 

. B.l  Geometry 

The  sensor  is  carried  in  a satellite  in  some  orbit  and  is 
allowed  to  scan  either  side  of  the  orbital  plane.  The  scan  pattern 
lies  on  a surface  of  a cone  whose  axis  goes  through  the  center  of 
the  earth  and  the  satellite.  This  geometry  gives  a constant  range 
to  the  ground  clutter  for  a specific  half  cone  angle.  The  geometry 
is  shown  in  Fig.  B-1.  The  velocity  vector  of  the  satellite  is 
normal  to  the  plane  of  Fig.  B-1. 

The  target  will  be  at  some  altitude  which  is  different 
than  the  altitude  of  the  satellite.  Figure  B-2  shows  the  plane 
AOB  and  defines  the  terminology  of  various  angles  referred  to  in 
this  note.  All  angles  are  measured  to  the  center  of  the  antenna 
beamwidth. 
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The  ground  clutter  returns  will  be  spread  over  a range 
interval  defined  by  the  intersection  of  the  antenna  beanvwidth  and 
the  earth.  Figure  B-3  is  an  amplitude  versus  time  display  (A-scope) 
showing  the  problem  analyzed. 

The  target  return  may  or  may  not  fall  within  the  ground 
clutter  return.  Its  position  is  a function  of  antenna  beamwidth, 
depression  angle,  satellite  altitude,  and  target  altitude.  Clearly, 
if  the  target  is  on  the  ground  it  will  be  imbedded  in  the  clutter 
and  other  techniques  must  be  used  to  separate  the  target  and  the 
clutter.  However,  when  the  target  is  above  the  surface  of  the 
earth,  there  is  a combination  of  the  above  mentioned  parameters 
which  result  in  the  target  being  at  a range  closer  than  the 
minimum  range  of  the  clutter. 

B.2  Mathematical  Development 

Figure  B-4  defines  the  terminology  in  the  AOB  plane  shown 

in  Fig.  B-1.  The  two  quantities  that  must  be  solved  for  are  the 

maximum  target  range,  R.  , and  the  minimum  range  to  the  ground 

^max 

clutter.  These  cwo  occur  at  opposite  edges  of  the  antenna  beam. 

The  minimum  range  to  the  clutter  is  given  by  AC.  From  the  law  of 
sines  for  triangle  AOC 
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Fig.  B-3.  A-scope  display. 
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where 


Y- 


7r/2  - + BW/2) 


BW 


depression  angle 


3 dB  radar  beamwidth 


a- 


earth  central  angle 


. -1 


= Ti-sin 


f/%  . 1 

^ n -^1 

L e J 


(2) 


In  a similar  manner  the  maximum  range  to  the  target  is 


given  by 


(sin  a, 

SIP. 


\ = R, 

max 


(3) 


where 


Y2  = Tr/2  - (ip  - BW/2) 


. -1 


* hx  1 

A -J 


^2  " ^"sxn 


(4) 


The  range  separation  between  clutter  and  target  is  given 


by 
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AR  = R ~ \ 

^min  max 

[sin  _ sin  02 
sin  y^  sin 

Positive  solutions  of  equation  6 indicate  when  the  target 
is  at  a range  less  than  the  clutter  and  range  gating  can  be  used 
to  eliminate  the  clutter.  Negative  solutions  of  equation  6 indicate 
regions  where-  Doppler  filtering  must  be  employed  to  separate  the 
target  from  the  clutter.  Equation  6 has  been  evaluated  as  a 
function  of  the  four  variables;  satellite  altitude,  target  altitude, 
depression  angle  and  3 dB  antenna  beamwidth.  The  results  are  pre- 
sented in  the  next  section. 

A satellite-borne  search  system,  to  be  cost  effective, 
should  cover  as  large  a portion  of  the  earth's  surface  as  possible, 
otherwise  the  system  will  require  a large  number  of  satellites. 

The  ground  swath  width  in  this  system  depends  on  the  satellite 
altitude,  depression  angle  and  scan  angle.  Consider  the  plane 
defined  by  the  intersection  of  the  cone  scan  pattern  and  the 
earth  as  shown  in  Fig.  B-5  for  which  the  following  definitions 
apply. 

EF  is  a great  circle  arc  whose  length  is  2R^a 
a is  the  earth  central  angle  defined  in  Fig.  B-1 
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GH  is  a great  circle  arc  whoso  length  is  the 
desired  ground  swath  width 

^KJG  = 6 the  required  half  scan  angle 

Given  the  above,  the  components  of  triangle  KJG  can  be 
determined  to  be 


KG  = 


e 


(7) 


JG  = R sin  a 
e 


(8) 


The  required  half  scan  angle,  6 , to  achieve  a particular 
swath  width  can  be  determined  from 


sin  a 


Equation  9 is  used  in  the  next  section  to  determine  the 
swath  width  limits. 

B.  3 Discussion  of  Results 

In  this  section  the  results  of  the  analysis  are  presented 
in  graphical  form.  The  range  separation  between  target  and  clutter 
is  plotted  as  a function  of  depression  angle.  Satellite  altitude, 
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3 dB  radar  beamwidth  and  target  altitude  are  parameterized. 

Figures  B-6  through  B-8  show  the  range  separation  versus 
depression  angle  for  a target  at  10  km  altitude  and  beamwidths  of 
0.1,  0.2  and  0.3  degree  respectively.  Satellite  altitude  is 
parameterized  on  each  curve  from  1000  to  6000  km.  There  are  two 
limits  shown  on  each  curve.  These  are  the  earth's  limb  which 
occurs  at  small  depression  angles  and  an  arbitrary  1000  nautical 
mile  swath  width  limit.  Operation  on  the  dashed  portion  of  the 
curves  is  possible  at  smaller  swath  widths. 

Figures  B-9  and  B-11  show  similar  results  for  targets  at 
8 km  altitude  and  Figs.  B-12  through  B-14  show  the  results  for 
targets  at  6 km  altitude. 

These  figures  indicate  that  range  gating  may  be  used  for 
a large  range  of  satellite  altitudes  against  targets  at  altitudes 
to  6 km  and  below.  As  the  3 dB  radar  beamwidth  is  increased,  the 
region  where  this  technique  is  applicable  diminishes.  For  example, 
at  satellite  altitudes  above  1000  km  and  a 3 dB  radar  beamwidth  of 
0.3  degree.  targets  below  6 km  altitude  will  be  at  ranges  comparable 
with  the  ground  clutter  if  a 1000  n.mi.  wide  ground  swath  must  be 
covered. 

B.4  Summary 

An  analysis  has  been  performed  to  show  the  feasibility  of 
using  range  gating  as  a technique  to  eliminate  main  lobe  ground 
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Range  separation  vs.  depression  angle. 
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Fig.  B-7.  Range  separation  vs.  depression  angle. 
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Fig.  B-8.  Range  separation  vs.  depression  angle. 
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Fig.  B-9.  Range  separation  vs.  depression  angle. 
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Fig,  B-12 . Range  separation  vs.  depression  angle. 
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clutter  from  the  radar  return  in  a satellite-borne  radar  search- 
ing an  altitude  regime  above  the  surface  of  the  earth.  A wide 
range  of  satellite  altitudes,  target  altitudes,  sensor  depression 
angle  and  beamwidths  were  investigated.  This  technique  appears 
to  be  feasible  for  target  altitudes  above  6 km  and  3 dB  radar 
beamwidths  smaller  0.3  degree.  Operation  against  lower  altitude 
targets  is  possible  at  the  expense  of  covering  a smaller  ground 
swath . 
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